Topology and the Cosmic Microwave 

Background 



Janna Levin 



DAMTP, CMS, Cambridge University, 
Wilberforce Rd, Cambridge CBS OWA, U.K. 
j. levin@damtp. cam. ac. uk 
and 

Theoretical Physics Croup, Imperial College, London 



Abstract 

Nature abhors an infinity. The limits of general relativity are often signaled by 
infinities: infinite curvature as in the center of a black hole, the infinite energy of 
the singular big bang. We might be inclined to add an infinite universe to the list of 
intolerable infinities. Theories that move beyond general relativity naturally treat space 
as finite. In this review we discuss the mathematics of finite spaces and our aspirations 
to observe the finite extent of the universe in the cosmic background radiation. 
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1 Introduction 



Is the universe infinite? The universe had a birth and will have a death, whether from 
old age or a more catastrophic demise through a big crunch. With the advent of general 
relativity, cosmology has accepted space as an evolving geometry. Although the universe is 
believed to begin and end a dynamical life, we continue to ascribe a property to the universe 
that we would never assign to any other natural object, and that is the property of being 
spatially infinite. 

The universe is home to a plethora of structures all of which must be finite. It would be 
untenable to suggest that any physical structure be infinite, yet cosmologists often cavalierly 
assume space itself is infinite. Yet an infinite universe means an infinite number of stars 
and galaxies and people. By the law of probability an infinite universe accommodates an 
infinite number of events each happening an infinite number of times |2|. Somewhere else 
in the cosmos, you are there. In fact there are an infinite number of you littering space. An 
infinite number of us living the same lives and an infinite number living slightly different 
lives and all range of variants. While this concept of infinity is manageable to some, it 
reeks of something pathological to others. Is the universe a surreal object, unlike any of its 
inhabitants, or is it real and physical and limited and finite? 

If the universe is finite, then light will wrap around the finite space decorating the sky 
with ghost images, sometimes called clone images. In principle we should be able to see 
copies of ourselves at different ages in different directions. A kaleidescopic version of 'This 
is your life'. In practice, the universe does not appear to be so small. The universe, if finite, 
must at least be big. Big enough to house clusters of galaxies and perhaps even bigger than 
the entire light travel time since the big bang. Therein lies our question. If the universe is 
finite, how small is it, and how do we measure its shape? 

The question of the infinite extent of space is often answered dynamically by Einstein's 
theory. If the amount of matter and energy is underdense then the universe will expand 
forever, if overdense then the universe will ultimately recoUapse to a big crunch, and if 
critical then the universe will expand forever. To each of these dynamical possibilities 
corresponds a geometry of spacetime: the underdense cosmos is negatively curved and 
infinite, the overdense cosmos is positively curved and finite and the critical cosmos is flat 
and infinite. 

Each of these three geometries can support topologies with finite volumes without al- 
tering the dynamics or the curvature. Multiconnected topologies are often overlooked in 
favor of the simply connected possibilities. This oversight is encouraged by a limitation of 
Einstein's theory which elevates gravity to a theory of geometry but does not provide a 
theory of topology. If unification of gravity with the other forces ever succeeds, topology 
will inevitably be integrated as a predictive feature of any cosmology. The topology of 
additional compact dimensions is already understood to be a crucial feature in superstring 
theories and, in general, discretization of space into topologically finite bundles is a natu- 
ral consequence of quantization |Q. No one is good enough at quantum gravity or string 
theory to make a prediction for the geometry and shape of the whole universe. Still the 
requirement of extra topologically compact dimensions has become a commonplace notion 
in fundamental physics even if it seems esoteric in astronomy. The first suggestion of com- 
pact extra dimensions began with Kaluza back in Einstein's day and have been dressed up 
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in the modern guise of superstring theories since. Even if we can never see the extent of 
the large three dimensions, we may be able to look for small extra dimensions to find clues 
about the geometry of the cosmos. 

In the absence of a predictive theory we can still explore the mathematical possibilities 
for a finite universe. This review covers topologies consistent with the standard cosmological 
geometries and our aspirations to observe the global extent of space in patterns in the oldest 
historical relic accessible, the cosmic microwave background (CMB). Methods to determine 
topology from the CMB range from the geometric such as the search for circles in the 
sky 1 12, ^and topological pattern formation |14, 15, 16 1 to the brute force statistical 
methods such as the method of images [|l^, These methods will be explained in turn. 
This review is intended both for mathematicians interested in cosmology and cosmologists 
interested in the mathematics of topology. Not much will be taken for granted and we will 
take the time to explain even the standard cosmology. 

Topology can be observed with the CMB or with the distribution of collapsed objects 
such as quasars or clusters of galaxies. There are reviews on galactic methods ||3|, Q. To 
keep the scope of this article manageable we will limit our discussion to searches for artifacts 
of topology in the CMB and defer a discussion of galactic methods to the reviews in Ref. 
1^, P]. (see also Refs. [P, 0].) Also, much of the history on topology can be found in 
an earlier review |||]. Other accessible reviews can be found in Refs. |^, Collections of 
papers for two international workshops can be found in |1C] and [^]. The observational 
section in this review will focus on recent methods only, starting where Ref. |||] left off. 
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2 Topology and Geometry 



Cosmic topology aims to deduce the global shape of the universe by experimentally observ- 
ing a pattern in the distribution of astronomical objects. To do so, it is helpful to first 
understand the possible topologies from a purely mathematical perspective. The isotropy 
of the cosmic microwave background radiation implies that the curvature of the observable 
universe is very nearly constant, so in the present article we consider only manifolds of con- 
stant curvature. The homogeneity of the observable universe does not, of course, exclude 
the possibility that the curvature of space varies enormously on a global scale beyond our 
view, but because that hypothesis is presently untestable we do not pursue it here. By 
the same reasoning, we are most interested in multiconnected spaces that are small enough 
to witness, although theories beyond general relativity, such as string theory, may help us 
to push beyond those limitations (see §|^. There is a previous Physics Reports on cos- 
mic topology which provides an excellent overview of the mathematical principles The 
reader is encouraged to consult that article as well as Refs. [19, 120]. We will not repeat 



these detailed reviews but do survey the same topological methods for completeness and 
extend the discussion to include some additional modern methods. 

The global shape of any space, including the ultimate outerspace, is characterized by 
a geometry and a topology. The term geometry, as used by mathematicians, describes the 
local curves while topology describes global features which are unaltered by smooth defor- 
mations. A continuous transformation is known as a homeomorphism, a smooth continuous 
and invertible map which deforms one manifold into another without cutting or tearing. 
General relativity is not invariant under homeomorphisms but is invariant under diffeomor- 
phisms, which amount to a change of coordinates. It is this covariance of classical gravity 
which underlies the principles of general relativity. Covariance ensures that all observers 
experience the same laws of physics regardless of the worldline along which they travel 
and therefore regardless of the coordinate system within which they interpret the world. 
Relativity does not invoke such a principle with respect to topology. For this reason we 
do not have a fundamental principle to guide us when contemplating the topology of the 
universe. We might hope that such a principle or at least a prediction might precipitate 
from a complete theory of gravity beyond classical. In the meantime, we consider any 
mathematically allowed topology, restricting ourselves to the constant curvature spaces as 
indicated by observational cosmology. 



2.1 Overview of Geometry 

A cosmological principle asserts that the Earth should not be special in location or perspec- 
tive. Cosmological observations on the largest scales indicate we live in a space with the 
symmetries of homogeneity, space is the same in all directions, and isotropy, space looks 
the same in all directions. Homogeneous and isotropic manifolds have constant curvature. 
The n-dimensional geometries of constant curvature are the everywhere spherical S", flat 
and hyperbolic H" geometries (fig. ||). One natural consequence of curved geometries 
worth remembering is that the sum of the interior angles of a triangle are greater than vr 
on S", exactly vr on E", and less than vr on H". 

In 2D, all Riemannian surfaces are homeomorphic to the three constant curvature ge- 
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Figure 1: A homogeneous isotropic space may be spherical (positive curvature), flat (zero 
curvature), or hyperbohc (negative curvature). 



ometries S^, E^, and H^. The infinite hyperbohc plane cannot be drawn properly in 3D 
and so is diflicult to visualize. It can be drawn schematically as a surface which everywhere 
has the curvature of a saddle. A nice way to treat is as a pseudosphere embedded in 
(2 + l)-Minkowski space. The pseudosphere has the equation — + + = —1. The full 
isometry group of the hyperboloid is PSL{2,H) = 5L(2,R)/Z2 with 5L(2,R) the special 
Lorentz group of real 2x2 matrices with unit determinant. The embedding in Minkowski 
space makes the symmetry group of more obvious. 

The sphere can be embedded in Euclidean 3-space as a surface with radius + + 

= 1. The isometry group for the sphere is 0(3), all orthogonal 3x3 matrices with the 
absolute value of the determinant equal to one. 

The plane E^ is the infinite surface (x, y) with the Minkowski time coordinate z fixed 
at unity (see fig. |l^). The Euclidean plane has the full Galilean group of isometrics: trans- 
lations, refiections, rotations and glide refiections. A glide refiection involves a translation 
with a refiection along the line parallel to the translation. 

In 3D, all manifolds are not homeomorphic to the constant curvature manifolds. Instead 
there are 8 homogeneous geometries some of which are anisotropic. These 8 geometries 



were classified by Thurston in the mathematical literature [21, Cosmologists are more 



familiar with the equivalent classification of Bianchi into 8 homogeneous manifolds |23|. 
Out of respect for the cosmological principle we will consider the fully homogeneous and 
isotropic spaces of constant curvature S^,E^, and H^. Similar to 2D, can be embedded 
as a pseudosphere in (3 + l)-Minkowski space, as the sphere and E^ as the plane at 
fixed Minkowski time. This isometry group of S'^ is SO (A). The isometry group of E'^ is 
the product of translations and the special orthogonal 3x3 matrices. The isometry group 
of H3 is PSL{2,C) = SL{2,C)/Z2. 



2.2 Overview of Topology 

Topology is the pursuit of equivalence classes of spaces; that is, the classification of homeo- 
morphic manifolds. To a topologist, all surfaces with one handle for instance are equivalent. 
Thus the doughnut and the coffee cup, in addition to providing a dubious breakfast, pro- 
vide the quintessential example of equivalent topologies, despite their obviously different 
geometries (i.e. curves). 

The doughnut and coffee cup are both manifestations of the 2-dimensional torus (fig. |2|), 
a prototypical multiply connected space. The torus has a finite area, yet has no boundary. 
When an inhabitant of the torus looks forward (fig. P), her line of sight wraps around the 
space and she sees herself from behind. She has the illusion of seeing another copy of herself 
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Figure 2: The flat 2-dimensional torus is defined by abstractly gluing opposite edges of a 
square. 




Figure 3: An inhabitant of the 2-torus looks forward and sees herself from behind. 

directly in front of her (fig. ^). Indeed she also sees herself when she looks to the side, or 
along a 45 degree line, or along any line of rational slope. She thus has the illusion of seeing 
infinitely many copies of herself (fig. ^). In other words, even though her universe is flnite, 
she has the illusion of living in an infinite universe containing an infinite lattice of repeating 
objects. 

As described below, all multiconnected surfaces can be built by gluing the edges of a 
fundamental polygon. The torus is built by identifying opposite edges of a parallelogram. 
If we begin with a fiat rectangle, we can bend this fiat sheet into three-dimensions to glue 
the top to the bottom and the left edge to the right. In doing so we have made a torus 
of revolution as in fig. |^ which is neither flat nor constant curvature. The curvature of the 
torus clearly varies over the surface. This bending of the torus is an artifact of embedding 
the surface in 3-dimensions. A truly fiat torus, better known as T'^ = x S^, is content to 
live in 2D with no such bending and projecting into 3D. The flat torus is akin to the video 
game with periodic boundary conditions where the left edge is identified with the right edge 
so that a flat explorer could stick their hand out the right edge only to have it appear poking 
out of the left edge. The edge in fact is as meaningless to the inhabitant of a flat torus as it 
is to an inhabitant of the torus of revolution. is truly flat and so has a different geometry 
from the torus of revolution although they are topologically equivalent, being characterized 
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Figure 4: The inhabitant of the 2-torus has the illusion of seeing a copy of herself. 
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Figure 5: Indeed, the inhabitant of the 2-torus sees a lattice of images of herself. 
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Figure 6: The torus of revolution has the same topology as the flat torus but different 
curvature. 




Figure 7: Left: The Klein bottle is similar to the torus, only the square's top and bottom 
edges are each glued to themselves, not to each other. The lattice of images in the Klein 
bottle is different from that in the torus because it contains glide reflections as well as 
translations. 

by the one handle. The genus g of the manifold is a topological invariant which counts the 
number of handles and holes. The torus has genus g = 1- 

The 2-dimensional Klein bottle (fig. 0) is similar to the torus, only now the top and 
bottom of the rectangle are glued not to each other, but each to itself with a shift of half 
a unit. The traditional way to make a Klein bottle is to start with the rectangle, glue top 
to bottom but glue the left to the right after a rotation through vr. In the present article 
we use a different (but equivalent) construction (namely gluing the top and bottom each to 
itself with a shift of half a unit) because it will simplify the construction of flat 3-manifolds 
in section ^.2| . As in the torus an inhabitant has the illusion of seeing an infinite lattice 
of objects, only now the structure of the lattice is different: it contains glide reflections as 
well as translations. The truly flat Klein bottle is topologically equivalent to the bottle 
immersed in 3D and pictured in fig. |8[ The embedding again leads to a curved surface and 
in this case to a self-intersection of the surface. The flat Klein bottle, content to live in 2D, 
has no such curvature or self-intersection. 

These examples generalize readily to three dimensions. A cube (fig. ^ with opposite 
faces glued is a 3-dimensional torus, or 3-torus. Its inhabitant has the illusion of seeing 
an infinite 3-dimensional lattice of images of himself and of every other other object in the 
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Figure 8: A surface with the topology of the Klein bottle. 




Figure 9: 3-dimensional lattice of cubes. 

space. Varying the gluings of the faces varies the structure of the lattice of images. Cosmic 
topology aims to observe these distinctive lattice- like patterns. 

2.3 Tools of Topology 

In this preliminary discussion several important concepts have already begun to emerge. 
The first is the fundamental domain. A fundamental domain is a polygon or polyhedron 
from which a manifold may be constructed. For example, fig. ^ shows a square fundamental 
domain for a 2-torus, and fig. |^ shows a cube fundamental domain for a 3-torus. A Dirichlet 
domain is a special type of fundamental domain. To construct a Dirichlet domain, pick an 
arbitrary point in the manifold of interest to serve as a basepoint. Start inflating a balloon 
with center at the basepoint (fig. p!o| ) and let the balloon expand uniformly. Eventually 
different parts of the balloon will bump into each other. When this happens, let them press 
flat against each other, forming a flat (totally geodesic) boundary wall. Eventually the 
balloon will fill the whole manifold, at which point it will have the form of a polygon in two 
dimensions or polyhedron in three dimensions whose faces are the aforementioned boundary 
walls. This polyhedron is the Dirichlet domain. Gluing corresponding faces recovers the 
original manifold. Notice that the fundamental rectangular domain for the Klein bottle 
when glued after a shift of half a unit is not Dirichlet, but when glued after a rotation 
through TT is Dirichlet. In some manifolds the Dirichlet domain depends on the choice of 
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Figure 10: Start inflating a balloon, with its center point fixed, and let it expand until it 
fills the entire space. The resulting polygon (in two dimensions) or polyhedron (in three 
dimensions) is called a Dirichlet domain. In some manifolds the shape of the Dirichlet 
domain depends on the choice of the center point; in other manifolds it does not. 



basepoint while in others it does not - more on this later. 

By identifying the edges of the fundamental domain, a multiconnected space is con- 
structed. Figures H and ^ show how an inhabitant of a finite multiply connected space may 
have the illusion of living in an infinite simply connected space. This apparent space is 
called the universal cover and can be thought of as the simply connected manifold with the 
constant curvature geometry. The fundamental domain is cut from this cloth and glued 
together while preserving the geometry of the universal cover. So the flat torus preserves 
the geometry of the universal cover while the torus of revolution does not. 

The group of covering transformations is the group of symmetries of the universal cover. 
That is, it is the group of motions of the universal cover that take images of a given object 
to other images of the same object. For example, in the case of a torus the group of covering 
transformations is the obvious lattice of translations. In the case of the Klein bottle the 
group of covering transformations contains glide reflections as well as translations. 

In two-dimensions all topologies have been completely classified. We give that clas- 



sification in §3.1. However, in three-dimensions the classification of hyperbolic manifolds 
remains incomplete although huge advances were made last century. There is a collection of 
topological invariants which can be used to recognize the equivalence classes of spaces. We 
have already encountered the genus of the manifold which counts handles. A group struc- 
ture can be used to capture handles and thereby detect the connectedness of the space. 
The most important of these is the loop group also known as the fundamental group or 
the first homotopy group. All loops which can be smoothly deformed into each other are 
called homotopic. A loop drawn on the surface of a sphere can be contracted to a point. In 
fact all loops drawn on the surface of a sphere can be contracted a point. Since all loops 
are homotopic, the sphere is simply connected and in fact does not have any handles. If 



however we draw a loop around the torus of revolution as in fig. 11, there is no way to 



smoothly contract that loop to a point, similarly for a loop drawn around the circumference 



of the toroid also shown in fig. 11 . The torus is therefore multiply connected. The group 
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Figure 11: Two homotopically distinct loops around the torus of revolution. 



of loops is the fundamental group 7ri(7W) and is a topological invariant. In more than 2D, 
the fundamental group is not sufficient to determine the connectedness of the manifold. 
Nonetheless, the fundamental group is a powerful tool. Poincare has conjectured that if the 
fundamental group is trivial, then the n-manifold is topologically equivalent to S"". 

A closely related structure is the holonomy group. The holonomy group can be thought 
of as the set of instructions for identifying the faces of the fundamental domain. The 
holonomies act discretely and without fixed point. The holonomy group is a subset of the 
full isometry group of the covering space. The full symmetries of the universal cover are 
broken by the identifications and it is customary to represent the compact manifold as a 
quotient space G/T where G is the isometry group of the universal cover and T is the holon- 
omy group. The holonomy group is of extreme utility in the game of cosmic topology. The 
group provides the set of rules for distributing images throughout the universe and thereby 
provides the mathematical means by which to determine the geometric distribution of as- 
tronomical images in the cosmos. The group is synonymous with the boundary conditions 
for the manifold and as we will see this is critical for building a prediction for the spectrum 
of the cosmic microwave background. 

There can be many different representations for the holonomy group. There is often a 
simplest presentation {71. ..7n;?'} where 7^ S T are elements of the group. Implicitly r = 1 
following the semicolon and is a list of all of the relations among the group elements. The 
most natural presentation for cosmology has a group element, and so a letter, for every 
identification rule. Such a presentation has many letters and many relations. 

The holonomies can map the end-point of an orbit to the start of that orbit so that it 
becomes periodic in the compact space. In other words, a periodic orbit has Xg^d — 7''^start 
where 7 can be a composite word 7 = n^Tfei- Each word has a corresponding periodic 
orbit. 

As an explicit example we return to the canonical torus. We can make use of the 
embedding of in a (2 -|- l)-dimensional space. Specifically, the 2-dimensional coordinate 
is replaced with the (2 -|- l)-dimensional coordinate 



(2.1) 




where tm is fixed at unity as in fig. 12. In this coordinate system the holonomy group for 
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Figure 12: The embedding of into (2 + l)-Minkowski space. The manifold appears as 
an infinite sheet fixed at tm = 1- 




Fi gure 13i A periodic orbit on the torus which corresponds to sjen^j 
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The boundary condition is then x'^ 
right edge results in 
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Tl^x^ so that identification of the left edge to the 



(2.3) 





The holonomy group for the torus in this presentation is {Tx.,Ty;TxTyT^^Ty^^. The relation 
TxTyT~^Ty^ = 1 is simply the commutivity of the generators and prunes the number of 
unique periodic orbits and therefore the number of clone or ghost images. Because the 
generators commute, the periodic orbits (homotopic loops) can be counted symbolically in 
terms of the number of windings taken around x and around y, {nix^rny). The number of 
loops and therefore the number of clone images grows as a polynomial with length. For 
example, consider the periodic orbit of fig. 13 which has Xend = TyT'^Xst&rt and so winding 
number (m^, = 2,my = 1). 

The non-orientable Klein bottle first twists the ^-faces through tt before identification. 
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The generators of F for the twisted space are Tx,R{'K)Ty with R{0) the rotation matrix: 



R{e) 



( 



1 

cos 9 sin 6 
— sin 9 cos 9 



) 



(2.4) 



All of the multiconnected, flat topologies can be built out of a combination of rotations and 
these translations. 

In hyperbolic space, the generators cannot easily be written in functional form but 
can nonetheless be found with numerical entries in the matrix. The program SnapPea 
provides the generators explicitly in this embedded coordinate system for a census of known 
hyperbolic manifolds and orbifolds (Orbifolds can have singular points.) The number 
of periodic orbits and therefore of ghost images grows exponentially in a compact hyperbolic 
manifold. The exponential proliferation of images is a manifestation of the chaotic flows 
supported on compact hyperbolic spaces and discussed in § |7.3| . 

2.4 Tessellations 

The modest tools we have introduced go a long way in constructing cosmological models. 
To summarize we will characterize a manifold by: 

• The fundamental domain - the shape of space around a given basepoint. 

• The universal cover - the simply connected manifold with the same constant cur- 
vature geometry. 

• The holonomy group - face pairing isometrics, i.e., the rules for identifying the faces 
of the fundamental domain. 

Armed with these tools, we can begin to visualize life in a finite universe. An important 
visualization technique comes from a tessellation of the universal cover. Consider the torus 
again. The flow of an observer or of light can be followed through the compact flat space 
by conscientiously respecting the boundary conditions and drawing the motion within the 
fundamental polygon as in fig. The Dirichlet domain is the set of points 



In T^, light travels along straight lines as do all inertial observers. Another faithful represen- 
tation of T'^ which does not distort the everywhere flat nature of the geometry is provided 
by a tiling picture. Beginning with the fundamental polygon, the periodic boundary con- 
dition which identifies the left edge of the rectangle to the right edge for instance can be 
represented by gluing an identical copy of the fundamental polygon left edge to right edge. 
This amounts to moving an entire copy of the fundamental domain with Tx{FD). Similarly 
identical copies can be glued top to bottom Ty[FD) ad infinitum until the entire infinite 
flat plane is tiled with identical copies of the rectangle. In general, the generators of the 
holonomy group can be thought of as defining an alphabet of 2n letters, for n generators and 
their n inverses. The alphabet will vary with different representations, as will the grammar. 
Usually, though not always, the longer the word built out of the letters 7, £ T, the more 



(2.5) 



such that 




(2.6) 
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Figure 14: The compact hypertorus can be represented as an identified parallelepiped. Al- 
ternatively, the compact topology can be represented by tiling flat space with identical copies 
of the fundamental parallelepiped. In the tiling picture above only the {x,y) directions are 
shown. A particular periodic orbit is drawn which corresponds to Xend = JyT^Xstart- 

distant the corresponding image or tiling in the tessellation picture. Figures ^, ^, and 
demonstrate tessellations. 

These copies are truly identical. If someone stands in the center of one rectangle he will 
see images of himself standing at the center of every other rectangle. If he moves to the left 
he will see every copy of himself move left. Because of the finite light travel time, it will take 
longer to see the more distant images movings. The light from these more distant images 
takes a longer path around the compact space before reaching the observer, the source of 
the very image. This creates the pattern of images already described so that the universe 
is an intricate hall of mirrors. 

Notice that the tiling completely fills the plane without overlaps or gaps. We can also 
imagine tiling the plane with a perfect square or even a parallelogram without overlaps or 
gaps. A fiat torus can therefore be made by identifying the edges of these other polygons. 
All have the same topology and the same curvature but some have different metrical features 
in the sense that the rectangle is truly longer than it is wide. An observer could line rulers 
up in space and unambiguously measure that their space was longer than it was wide while 
an observer in the perfect square would know that their space was equally long and wide. 
The torus is thus globally anisotropic even though the geometry of the metric is still locally 
isotropic. The torus is still globally homogeneous. The center could be moved arbitrarily 
and the fundamental polyhedron redrawn around that center to give an equally valid tiling. 
Most topologies globally break both homogeneity and isotropy. 

The holonomies provide the list of rules for identifying the edges of the polygon. So the 
Klein bottle, while it can begin with the same fundamental polygon as the torus and has the 
same universal cover, it has a different set of rules for identifying the edges and therefore a 
different holonomy group. The tiling will be correspondingly different with identical copies 
glued top edge to bottom edge but left edge to right only after rotating the tile by vr. So 
the tessellation creates the refiected and translated images specific to the Klein bottle (fig. 

0). 

The fiat plane can also be tiled by hexagons (fig. |lO|) as many a bathroom fioor demon- 
strates. The hexagon therefore also reproduces a fiat torus. However the general parallelo- 
gram and hexagon exhaust the possible convex polygons which can tessellate fiat space. If 
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Figure 15: 



we tried to fill the plane with octagons, we would find that they overlapped at a vertex and 
could not smoothly fill the plane, as in fig. 15. If we glue the octagons together without 
confining them to the plane, they will create a curling surface that tries to have a hyperbolic 
structure. In fact, the 2D hyperbolic plane H'^ can be tiled by regular octagons. Recall 
that the angles of a triangle sum to less than tt in hyperbolic geometry. Since the interior 
angles of the polygon narrow on the hyperbolic plane, they can be drawn just the right 
size relative to the curvature scale so that precisely the right number fit around a vertex 
to fill the negatively curved plane without gaps or overlaps. Notice also, that there is no 
scale in flat space and there is no intrinsic meaning to the area of the compact torus. It 
can be made arbitrarily large or small relative to any unit of measure. This is not the case 
when the surface is endowed with a curved metric. The area is then precisely determined 
by the topology. The area of the octagon has to be precise relative to the curvature scale 
to properly tessellate the universal cover. 

The topology created by the compact octagon is a double torus T2 where the 2 refers 
to the number of handles or genus, g, of the manifold. Hyperbolic T2 is topologically 
equivalent to the double doughnut of fig. 16 which lives in although again T2 has a 
different geometry, being flat, while the double doughnut is clearly curved. 

A connected sum of tori can be constructed to make Tg, a flat 2D compact surface with 
g handles. In order to make a g handled object, a polygon with 4:g edges is needed. The 
polygons have to be drawn larger relative to the curvature scale the more edges they have 
in order that the interior angles be thin enough to fit together around a vertex and tile 
H^. Consequently a relationship emerges between the area of the surface and the topology. 
The area of the hyperbolic surface can be related to the Poincare-Euler characteristic. 
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Figure 16: 



X = 2(1 — g) through the Gauss-Bonnet theorem 

A = -27Tx = M9-l)- (2.7) 

Although the area is a topological invariant in 2D, the lengths need not be fixed. For 
instance, the asymmetric octagon and the regular octagon have the same topology, namely 
that of the double doughnut, and the same negative curvature, but they have a different 
spectrum of periodic orbits. 

The tessellation of the universal cover is of particular importance when we consider 
three-dimensional spaces. We simply do not have four dimensions available to us in which 
to bend the three-dimensional volume and that visualization technique fails us. However, 
we can still visualize filling the three-dimensional volume with 3D tiles which preserve 
the local geometry of the manifold. Therefore we rely heavily on the universal cover, the 
fundamental polyhedron and the holonomy group when investigating the topologies of these 
compact spaces. 

2.5 Hyperbolic topologies and Volumes 

All 2-dimensional manifolds have been classified, as have flat and spherical 3-manifolds, but 
hyperbolic 3-manifolds are not yet fully understood. These continue to resist classification 
although enormous advances were made last century. An important realization was the 
rigid connection between metrical quantities and topological features in hyperbolic three- 
manifolds. Metrical quantities describe lengths and would ordinarily fall under the domain 
of geometry. Yet in three-dimensions the Mostow [p5| ] rigidity theorem ensures that once 
the topology is specified, then all metrical quantities are fixed for 3-D hyperbolic manifolds. 
This means that not only the volume but also the lengths of the shortest geodesies are 
immutable for a given topology. This is to be compared with 2D where a given topology 
can support an infinite number of metrics as already described. 

Although the topology fixes the compact hyperbolic volume, the volume does not 
uniquely specify the topology. In other words there can be different topologies with the same 
volume. Because of this powerful link between geometry and topology, compact hyperbolic 
manifolds can be ordered according to their volume. A remarkable result is that compact 
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hyperbolic manifolds form a countably infinite number of countably infinite sequences or- 
dered according to volume. A given sequence shows an accumulation of manifolds near a 
limiting volume set by a cusped manifold. A cusped manifold has finite volume but is not 
compact, having infinitely long cusped corners which taper infinitely thin. The cusps are 
topologically a 2-torus which is conformally shrunk to zero down the narrowing throat. The 
sequence of compact manifolds are built systematically through the process of Dehn surgery 
which is used in the next section to organize the spectrum of topologies. Dehn surgery is 
a formal process whereby a torus is drilled out of the cusped manifold and replaced with a 
solid torus. The surgery is taken along a periodic orbit of identified by the number of 
windings {nix, my) the orbit takes around the torus before closure. The program SnapPea 
names manifolds according to the number of windings so that the manifold m003(3,-l) is a 
compact space constructed by Dehn filling the orbifold m003 along the periodic orbit with 
winding numbers (3,-1). 

The space m003(3,-l) is also known as the Weeks space after the mathematician who 
discovered it and coincidentally has contributed to this article [^^. The Weeks space is of 
particular relevance to cosmology since it is the smallest known manifold with a volume 
V ~ 0.94 in curvature units. The Weeks space displaced the Thurston space m003(-2,3) 
which has volume V ~ 0.98. Although the Weeks space is the smallest known manifold it 
may or may not be the smallest possible. The current bound on the minimum volume for 
any compact hyperbolic manifold is about 0.3 in units of the curvature radius. Contrast 
this with the flat hypertorus which can be made infinitesimally small. Again, the rigidity 
of compact hyperbolic spaces is at work. 

Because of the rigidity of 3-D hyperbolic manifolds, we can also characterize a topology 
by certain volume measures. Very useful quantities in cosmological investigations are the 
in-radius r_ and the out-radius r_|_. The in-radius is the radius of the largest geodesic ball 
centered at the origin which can be inscribed within the Dirichlet domain. The out-radius 
is the radius of the smallest geodesic ball centered at the origin which can encompass the 
Dirichlet domain. The disparity between r_ and + can be taken as a rough indication of the 
global anisotropy introduced by the topological identification. The injectivity radius is the 
radius of the largest ball, centered at the given point, whose interior embeds the manifold. 
The Dirichlet domain given by default in SnapPea is centered at a local maximum of the 
injectivity radius. In this way it often provides the most symmetric form for the Dirichlet 
domain. Moving the observer away from this basepoint will change the appearance of 
the global shape of space observationally. A somewhat ambiguous but sometimes useful 
quantity is an estimate of the diameter of the manifold. Since the manifold is not isotropic 
this is not precisely defined. Nonetheless, one estimation is given by 



that is, it is the largest distance between any two points in the Dirichlet domain. Notice 
that it is not a geodesic distance so no periodic orbit lies along the diameter. A more well 
defined diameter is the maximum distance between two points in the manifold, where the 
distance is measured in the manifold itself, not in the Dirichlet domain: 



dM = supx,yeMd{x,y); 




diam = supx,y{ 




(length of path)). 



(2.9) 
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The periodic geodesies are quite important and are closely related both to the eigen- 
modes in the space and the distribution of ghost or clone images. The periodic geodesies 
can be found with the holonomy group. The eigenvectors of the holonomy group span a 
plane in the AD embedding which intersects the pseudosphere H^. The line of intersection 
occurs along the periodic orbit corresponding to that group element. A composite element 
can be thought of as a word formed from the fundamental alphabet of the presentation. 
The fact that the number of words grows exponentially with length means an exponential 
growth of long periodic orbits. The proliferation of periodic orbits, dense in phase space, are 
the canonical mark of chaotic dynamics. This leads us to remark that compact hyperbolic 
spaces are not just complicated mathematically but they are truly complex, supporting 
chaotic flows of geodesies. Since most cosmological investigations try to circumvent the 
chaotic nature of the motions, we mention chaos only occasionally as we go along. We 
reserve a final section for a brief discussion of this fascinating feature. A review of chaos on 
the pseudosphere can be found in Ref. [27|. 
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Figure 17: To add a handle to a surface, remove two disks and glue the remaining boundary 
circles to each other. 



3 Survey of Compact Manifolds 

Having described the topological principles at work and some examples, we can compile 
a list of known manifolds, their geometries and topologies. The classification is clearly 
presented in Refs. ||l|, 20|. We present a different mathematical approach for variety. 



The classification is organized by J. Weeks |28] in relation to four questions: 

1. What constant curvature manifolds exist? 

2. How flexible is each manifold? That is, to what extend can we change a manifold's 
shape without changing its topology or violating the constant curvature requirement? 

3. What symmetries does each manifold have? 

4. Does a Dirichlet domain for the manifold depend on the choice of basepoint? 

It turns out that the answers to these questions depend strongly on both the curvature 
and the dimension. The remainder of §^ answers these questions in each of the six possible 
cases: that is, in spherical, flat, or hyperbolic geometry in two dimensions (§^.l|) or three 
dimensions (§ |3.2[) . 



3.1 Two-dimensional manifolds 

3.1.1 Overview of two-dimensional manifolds 

Every closed surface may be obtained by starting with a sphere and adding handles (fig. 



17) or crosscaps (fig. 18). Each handle is constructed by removing two disks from the 
sphere and gluing the resulting boundary circles to each other (fig. [Tt] ), and each crosscap 
is constructed by removing one disk and gluing the resulting boundary circle to itself (fig. 

18) . Conway found a particularly simple and elegant proof of this classification; for an 
illustrated exposition see Ref. [p9| . 

The interaction between geometry and topology is a primary theme in our understanding 
of 2-manifolds. The promised interaction between geometry and topology is that each closed 
topological surface may be given a constant curvature geometry. Each closed 2-manifold 
admits a unique geometry. That is, a surface which admits a spherical geometry cannot also 
admit a flat geometry, and so on. For an elementary proof see Ref. |3^]. Infinite 2-manifolds, 
by contrast, may admit more than one geometry. For example, the Euclidean plane and the 
hyperbolic plane are topologically the same 2-manifold. We organize the surfaces according 
to the total number of disks that get removed during the construction (fig. [l9[) . 

In the case that no disks are removed, we have the sphere. 
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Figure 18: To add a cross cap to a surface, remove one disk and glue the remaining boundary 
circle to itself by identifying antipodal points. The construction cannot be physicahy carried 
out in without self-intersections, but abstractly the resulting surface has a seam that 
looks locally like the centerline of a Mobius strip. 

In the case that one disk is removed, we draw the sphere-minus-a-disk as a hemisphere, so 
that when the boundary circle gets glued to itself the spherical geometry continues smoothly 
across the seam from one side to the other, giving the resulting sphere- with-a-crosscap a 
uniform spherical geometry. 

In the case that two disks are removed we may have either a sphere- with-a-handle, 
better known as a torus, or a sphere-with-two-crosscaps, better known as a Klein bottle. 
We draw the sphere- minus- two-disks as a cylinder. The cylinder has an intrinsically flat 
geometry. You can construct one from a sheet of paper without stretching, and when you 
glue its top and bottom of the rectangle to each other (for the torus) or each to itself (for 
the Klein bottle) the flat geometry continues uniformly across the seams. Thus the torus 
and the Klein bottle each has a flat geometry. 

In the case that three disks are removed we have a sphere-with-three-crosscaps, con- 
structed from a three-way cylinder better known as a pair of pants. The pair of pants may 
in turn be constructed from two hyperbolic hexagons (fig. |20|) , giving it an intrinsically hy- 
perbolic geometry. When the pants' boundary circles (the cuffs) are glued to themselves to 
make the sphere-with-three-crosscaps, the hyperbolic geometry continues uniformly across 
the seams. Thus the sphere-with-three-crosscaps also gets a hyperbolic geometry. 

All remaining surfaces - those for which four or more disks are removed in the construc- 
tion - may be built from pairs of pants, so all get a hyperbolic geometry. In summary, we 
see that of the infinitely many closed 2-manifolds, two admit spherical geometry (the sphere 
and the sphere- with-crosscap) , two admit flat geometry (the torus and the Klein bottle), 
and all the rest admit hyperbolic geometry. 

The next three sections will examine the flat, spherical, and hyperbolic closed 2-manifolds 
in greater detail, in each case answering the four questions raised in § p.2| . 

3.1.2 Flat two-dimensional manifolds 

Manifolds: torus and Klein bottle. 
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Figure 19: Topological classification of closed 2-manifolds. Every closed 2-manifold may 
obtained from a sphere by adding either handles (left column) or cross caps (right column). 
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Figure 21: The torus has a two-parameter continuous family of symmetries, consisting of 
rotations and vertical translations. 



Figure 22: The Klein bottle has continuous rotational symmetry, but any attempt at vertical 
translation is blocked by the cross caps. 

Flexibility: The torus has three degrees of freedom in its construction: the height of the 
cylinder, the circumference of the cylinder, and the twist with which the top gets glued to 
the bottom. The Klein bottle has only two degrees of freedom in its construction: the height 
and circumference of the cylinder. Each end of the cylinder gets glued to itself, identifying 
diametrically opposite points, so there is no twist parameter. 

Symmetry: Continuous symmetries. The torus has a two-parameter continuous family 
of symmetries, consisting of rotations and vertical translations (fig. ^). This shows that the 
torus is globally homogeneous in the sense that any point may be taken to any other point by 
a global symmetry of the manifold. The Klein bottle has only a one-parameter continuous 
family of symmetries, consisting of rotations (fig. |2^ ). It is globally inhomogeneous. When 
diametrically opposite points on a boundary circle are glued to make a crosscap, the seam 
is a closed geodesic that looks locally like the centerline of a Mobius strip. Points lying on 
this centerline are fundamentally different from other points in the manifold. 

Discrete symmetries. As well as its continuous family symmetries, the torus has a 
discrete Z/2 symmetry given by a 180 degree rotation interchanging the cylinder's top and 
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Figure 23: The torus always has a discrete Z/2 symmetry given by a half turn, even if the 
cylinder's top and bottom circles are glued to each other with a twist. The Klein bottle has 
a discrete Z/2 + Z/2 symmetry, which interchanges and/or inverts the cross caps. 

bottom (fig. ^). The discrete and continuous symmetries may of course be composed. 
The best way to think of this is that the complete symmetry group of the torus consists 
of two disconnected components: symmetries that interchange the boundary circles and 
symmetries that do not. In special cases the torus may have additional discrete symmetries 
(fig. ^4|). The Klein bottle's discrete symmetry group is Z/2 + Z/2, generated by vertical 
and horizontal reflections of the cylinder. Vertical reflection interchanges the two crosscaps. 
Horizontal reflection takes each crosscap to itself, but reverses the direction of its centerline. 

Dirichlet domain: We have seen that the torus is globally homogeneous, so the shape of 
the Dirichlet domain does not depend on the choice of basepoint. On the other hand, the 
torus has three degrees of flexibility, and the Dirichlet domain's shape does depend on the 
shape of the torus. Typically the Dirichlet domain is an irregular hexagon, but in special 
cases it may be a square or a regular hexagon. The Klein bottle is globally inhomogeneous, 
so even for a fixed geometry of the Klein bottle, the Dirichlet domain's shape depends on 
the choice of basepoint (fig. p5| ). 

3.1.3 Spherical two-dimensional manifolds 

Manifolds: sphere and sphere-with-crosscap. 

Flexibility: The sphere and the sphere-with-crosscap have no flexibility whatsoever. 
Their geometry is completely determined by their topology. Here and throughout we make 
the assumption that spherical geometry refers to a sphere of radius one. Otherwise these 
two manifolds would have the flexibility of rescaling - they could be made larger or smaller. 

Symmetry: The sphere has a symmetry taking any point to any other point, with any 
desired rotation. These symmetries form a 3-parameter continuous family. It also has a 
single discrete symmetry, given by reflection through the origin. The sphere-with-crosscap 
has a 3-parameter continuous family of symmetries taking any point to any other point with 
any desired rotation and either reflected or not reflected. It has no discrete symmetries. 
Indeed the continuous family already takes a neighborhood of any point to a neighborhood 
of any other point in all possible ways, leaving no possibilities for discrete symmetries. 

Dirichlet domain: The sphere and the sphere-with-crosscap both have somewhat de- 
generate Dirichlet domains (try the inflating-the-balloon experiment in each), but these 
manifolds are globally homogeneous so the degenerate Dirichlet domains do not depend on 
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Cut the cylinder open along a line connecting a point on the bottom 
circle to the corresponding point on the top circle, taking into 
the twist in the gluing. The resulting fundamental domain is 
a parallelogram. 




Every parallelogram admit a half turn, for a Z/2 symmetry group. 
If the parallelogram happens to be a rectangle, it also admits 
reflections, for a D2 symmetry group. 

If the parallelogram happens to be a square, it admits a quarter 
turn, for a D4 symmetry group. 

Finally, if the parallelogram happens to be a rhombus with 
angle pi/3, the Dirichlet domain centered at a vertex is a regular 
hexagon, and the torus admits a one sixth turn, for a D6 symmetry group. 




Z/2 D2 D4 D6 



Figure 24: 




Figure 25: In a Klein bottle, the shape of the Dirichlet domain depends on the choice of 
basepoint. For example, a Dirichlet domain centered at the triple arrowhead is a hexagon, 
while a Dirichlet domain centered at a single arrowhead is a rectangle. 
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the choice of basepoint. 



3.1.4 Hyperbolic two-dimensional manifolds 



Manifolds: There are infinitely many closed hyperbolic 2-manifolds, illustrated in fig. 19. 

Flexibility: All closed hyperbolic 2-manifolds are flexible. The shape of each manifold 
is parameterized by the circumferences of the cuffs of the pairs of pants used to construct 
it, and by the twists with which distinct cuffs are glued to each other. A cuff that is glued 
to itself to make a crosscap has no twist parameter. 

Symmetry: Closed hyperbolic 2-manifolds never have continuous families of symmetries. 
The reason is that each seam (where cuffs are joined together) is the shortest loop in its 
neighborhood, and therefore cannot be slid away from itself. A seam cannot slide along itself 
either, because a pair of pants admits no continuous rotations. A generic closed hyperbolic 2- 
manifold also has no discrete symmetries, but nongeneric ones may have discrete symmetries 
if the cuff lengths and twists are chosen carefully. 

Dirichlet domain: The Dirichlet domain of a closed hyperbolic 2-manifold always de- 
pends on the choice of basepoint. For carefully chosen basepoints it may exhibit some of 
the manifold's symmetries. 

3.2 Three-dimensional manifolds 

3.2.1 Overview of three-dimensional manifolds 

Unlike two-manifolds, not all three-manifolds admit a constant curvature geometry. For 



a survey of the geometry of 3-manifolds, see Ref. |31]. Roughly speaking, any closed 3- 
manifold may be cut into pieces in a canonical way, and each piece admits a geometry that 
is homogeneous but not necessarily isotropic. For a very elementary exposition of homoge- 
neous anisotropic geometries, see Ref. [|^] For a more detailed but still readable exposition 
see 1 32]. The observed isotropy of the microwave background radiation implies that at least 



the observable portion of the universe is approximately homogeneous and isotropic, so cos- 
mologists restrict their attention to the three homogeneous isotropic geometries: spherical, 
flat and hyperbolic. However, as noted earlier, an approximately flat observable universe is 
also consistent with a very large universe of varying curvature. 

3.2.2 Flat three-dimensional manifolds 



Manifolds: Figure ^ illustrates the ten closed flat 3-manifolds. 

The first is the 3-torus introduced earlier - a cube with pairs of opposite faces glued. 

The second and third manifolds are variations on the 3-torus in which the front face 
is glued to the back face with a half (respectively quarter) turn, while the other two pairs 
of faces are glued as in the 3-torus. Other possibilities for gluing the faces of a cube - for 
example gluing all three pairs of opposite faces with half turns - fail to yield flat 3-manifolds 
because the cube's corners do not fit together properly in the resulting space. See Ref. [Q 
for a more complete explanation. 

The fourth and fifth manifolds are similar in spirit to the first three, only now the cross 
section is a hexagonal torus instead of a square torus. That is, the cross section is still 
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Figure 26: The ten closed flat 3-manifolds. 
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topologically a torus, but has a different shape (recall fig. [2^ ). The hexagonal cross section 
allows the front of the prism to be glued to the back with a one-sixth or one-third turn, 
producing two new flat 3-manifolds. The side faces are all glued to each other in pairs, 
straight across. 

The sixth manifold, the famous Hantzsche-Wendt manifold, is the most interesting of 
all. Unlike the preceding five manifolds, which we defined by constructing fundamental 
domains, the most natural way to define the Hantzsche-Wendt manifold is to start in the 
universal cover and define its group of covering transformations. Specifically, its group of 
covering transformations is the group generated by screw motions about a set of orthogonal 
but nonintersecting axes (indicated by heavy lines in fig. 24). Each screw motion consists 
of a half turn about an axis composed with a unit translation along that axis. Note that 
this group of covering transformations does not take a basic cube to all other cubes in the 
cubical tiling of 3-space. Rather, the images of a basic cube fill only half the cubes in 
the tiling, checkerboard style. Thus a complete fundamental domain would consist of two 
cubes, the basic cube and any one of its immediate neighbors; images of the basic cube 
would fill the black cubes in the 3-dimensional checkerboard, while images of its neighbor 
would fill the remaining white cubes. But we would really prefer a fundamental domain 
consisting of a single polyhedron. To get one we employ the balloon construction for a 
Dirichlet domain introduced in Section N.1.2. Let the Dirichlet domain's basepoint be the 
center of a basic cube. As the balloon expands its fills that basic cell entirely, and also 
fills one sixth of each of its six immediate neighbors. The resulting Dirichlet domain is a 
rhombic dodecahedron (fig. p^ ). The face gluings are given by the original screw motions 
along the axes. Note that this construction of the Hantzsche-Wendt manifold corrects an 
error, appearing elsewhere in the cosmological literature, which takes the three screw axes 
to be coincident. This erroneous construction leads to a cube with each pair of opposite 
faces glued with a half turn. The cube's corners are therefore identified in four groups of 
two corners each instead of a single group of eight corners. The resulting space has four 
singular points and is thus an orbifold instead of a manifold. 

The seventh through tenth manifolds are again defined by fundamental domains. They 
are similar to the first five manifolds, except now the cross section is a Klein bottle instead 
of a torus. That is, the fundamental domain's left and right faces are glued to each other 
as in a torus, but the top face (resp. bottom face) is glued to itself with a horizontal shift 
of half the width of the fundamental domain, taken modulo the full width, of course. This 
ensures that every cross section is a Klein bottle, which may be understood as a cylinder 
with one crosscap at the top and another at the bottom. The fundamental domain's front 
face, which is now a Klein bottle, may be glued to its back face in one of four ways: plain 
(seventh manifold), interchanging the crosscaps (eighth manifold), inverting the crosscaps 
(ninth manifold), or interchanging and inverting the crosscaps (tenth manifold). Unlike the 
first six manifolds, these four are all nonorientable. 

Flexibility: All ten closed flat 3-manifolds have at least two degrees of flexibility: the 
depth of each fundamental domain (fig. [2^ ) may be chosen independently of its width and 
height. For some of the manifolds, for example the one-sixth turn manifold, this is the only 
flexibility. Others have greater flexibility: for example the cross section of the half turn 
manifold may be any parallelogram, and the fundamental domain for the 3-torus may be 
any parallelepiped. 
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Symmetry: All closed flat 3-manifolds except the Hantzsche-Wendt manifold have a 
continuous one-parameter family of symmetries that pushes into the page in the fundamental 
domains in fig. These symmetries take each cross section to a cross section further along, 
with the cross sections that fall off the back of the stack reappearing cyclically at the front. 
Some of the manifolds have additional continuous families of symmetries as well. The 3- 
torus is the most symmetrical, with a 3-parameter family of symmetries taking any point to 
any other point. All ten manifolds have discrete symmetries as well, which vary from case 
to case and may depend on the manifold's exact shape as well as its topology. For example, 
a cubical 3-torus has a 3-fold symmetry defined by a one-third turn about one of the cube's 
long diagonals, but a 3-torus made from an arbitrary parallelepiped lacks this symmetry. 

Dirichlet domain: The 3-torus is exceptional in that its Dirichlet domain does not depend 
on the choice of basepoint. It is always the same, no matter what basepoint you choose, 
because there is a symmetry taking any point to any other point. For the nine remaining 
closed flat 3-manifolds, the Dirichlet domain always depends on the choice of basepoint. 

Observational status: All ten closed flat 3-manifolds are of course completely consistent 
with recent evidence of a flat observable universe. However, if the fundamental domain is 
larger than our horizon radius, the topology may, in practice, be unobservable. 



3.2.3 Spherical three-dimensional manifolds 

Manifolds: This section requires only a minimal understanding of the 3-sphere, or hyper- 
sphere, defined as the set of points one unit from the origin in Euclidean 4-space. Readers 



wanting to learn more about the 3-sphere may consult Ref. |3C] for an elementary exposition 
or Ref. ||2^ for a deeper understanding. 

Spherical manifolds fall into two broad categories. The manifolds in the first category 
have very different properties from those in the second. 

First category: 

The first category manifolds are most easily defined in terms of their groups of covering 
transformations, although it is not hard to work out their Dirichlet domains afterwards. 
Amazingly enough, their groups of covering transformations - which are symmetries of the 
3-sphere - come directly from finite groups of symmetries of the ordinary 2-sphere! In other 
words, each finite group of symmetries of the 2-sphere (that is, the cyclic, dihedral, tetrahe- 
dra, octahedral, and icosahedral groups) gives rise to the group of covering transformations 
of a spherical 3-manifold! 

The bridge from the 2-sphere to the 3-sphere is provided by the quaternions. The 
quaternions are a 4-dimensional generalization of the familiar complex numbers. But while 
the complex numbers have a single imaginary quantity i satisfying = —1, the quaternions 
have three imaginary quantities and k satisfying 

f = f^k^ = -i 
ij = k = -ji 
jk = i = —kj 

ki = j = -ik (3.10) 
The set of all quaternions a + bi + cj + dk (for o, b, c, d real) defines 4-dimensional Euclidean 
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space, and the set of all unit length quaternions, that is, all quaternions a + bi + cj + dk 
satisfying + + + (P = 1, defines the 3-sphere. 

Before proceeding, it is convenient to note that there is nothing special about the three 
quaternions i,j, and k. 

Lemma N.3.3.1 (quaternion change of basis). Let c be a 3x3 orthogonal matrix, and 
define 

3 ~ Cjii + Cjjj + Cjfjk 

k' = Ckii + Ckjj + Ckkk (3-11) 
then i', j' , and k' satisfy the usual quaternion relations 

i'j' = k' = -j'i' 
j'k' = i' = -k'j' 

k'i' = j' = -i'k' (3.12) 

Lemma N.3.3.1 says that an arbitrary purely imaginary quaternion bi + cj + dk may, by 
change of basis, be written as b'i' . If the purely imaginary quaternion bi + cj + dk has unit 
length, it may be written even more simply as i'. A not- necessarily-imaginary quaternion 
a + bi + cj + dk may be transformed to a' + b'i. If it has unit length it may be written as 
cos{9) + i' sm(9) for some 9. 

The unit length quaternions, which we continue to visualize as the 3-sphere, may act on 
themselves by conjugation or by left-multiplication. 

Lemma N.3.3.2 (conjugation by quaternions). Let q be an arbitrary unit length quater- 
nion. According to the preceding discussion, we may choose a basis {l,i' ,j' , k'} such that 
q = cos{9) + i' sin(6') for some 9. It is trivial to compute how q acts by conjugation on the 
basis {1, i',j', k'}. 

(cos(6l) + i' sin(6'))l(cos(6l) - i' sm{9)) = 1 

(cos(6') + i' sm{9))i' {cos{9) - i' sm{9)) = i' 

(cos(0) + i' sm{9))j'{cos{9) - i' s\n{9)) = j' cos(20) + k' sm{29) 

{cos{9) + i' sm{9))k'{cos{9) - i' s\n{9)) = -j' sui{29) + k' cos(20) (3.13) 

We see that conjugation always fixes 1 (the north pole), so all the action is in the equatorial 
2-sphere spanned by {i',j',k'}. The equatorial 2-sphere itself rotates about the i' axis 
through an angle of 29. 

Compare the preceding action by conjugation to the following action by left multiplica- 
tion. 

Lemma N.3.3.3 (left multiplication by quaternions). Let q be an arbitrary unit length 
quaternion. According to the preceding discussion, we may choose a basis {l,i' ,j' , k'} such 
that q = cos{9)+i' sin(^) for some t. It is trivial to compute how q acts by left multiplication 
on the basis {l,i' ,j' ,k'}. 

{cos{9) + i' sm{9))l = 1 cos{9) + i' sm{9) 
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(cos(^) + i'sin(^))i' = -1 sin(^) + i' cos(^) 
(cos(0) + i' sin(6))j' = j cos{0) + k' sm(0) 

(cos(^) +i'sin(0))A;' = -/sin(0) + A;'cos(0) (3.14) 

We see that left multiplication rotates 1 towards i' while simultaneously rotating j' towards 
k' . The result is a screw motion known as a Hopf flow. What is not obvious from the above 
computation is that the Hopf flow is completely homogeneous - it looks the same at all 
points of the 3-sphere. 

All the tools are now in place to convert finite groups of symmetries of the 2-sphere to 
finite groups of symmetries of the 3-sphere. The algorithm is as follows: 

Step ^1. Choose a finite group of symmetries of the 2-sphere. (The only such groups 
are the cyclic groups, the dihedral groups, the tetrahedra group, the octahedral group, and 
the icosahedral group. You might also expect a cubical and a dodecahedral group, but they 
coincide with the octahedral and icosahedral groups because the octahedron is dual to the 
cube and the dodecahedron is dual to the icosahedron.) 

Step #2. Write down the quaternions whose action by conjugation gives the chosen 
symmetries. Note that each symmetry corresponds to two quaternions q and —q, because 
q and —q act identically when conjugating, i.e. qrq"^ = {—q)r{—q)~^. 

Step ^3. Let the quaternions found in Step ^3 act by multiplication. This action 
defines the group of covering transformations. Note that although a quaternion's action by 
conjugation always has a pair of fixed points on the 2-sphere, its action by left multiplication 
on the 3-sphere is always fixed point free. 

The preceding algorithm is most easily understood in a concrete example. 

Example N.3.3.4. Consider the so-called Klein four group, the group of symmetries of 
the 2-sphere consisting of half turns about the i-, j-, and A;-axes. (Remember that we are in 
the 3-space spanned by {i, j. A;}, so the z-, j-, and A;-axes play the role of the traditional x-, 
y-, and z-axes.) By Lemma N.3.3.2 a half turn about the i-axis corresponds to conjugation 
by the quaternions ±i. That is, when 9 = 7r/2, the quaternion i = cos(7r/2) -|- isin(7r/2) 
acts by conjugation as a rotation about the i-axis through an angle 20 = ir. Similarly, the 
half turn about the j-axis corresponds to the quaternions ±j, and the half turn about the 
fc-axis corresponds to ib/c. The trivial symmetry corresponds to conjugation by ±1. Thus 
the complete set of quaternions is {ztl. zti, ztj, ztk}. If we now let these quaternions act by 
left multiplication on the 3-sphere, they give a group of covering transformations of order 8. 
(Exercise for the reader: compute the action of each of those eight quaternions on the basis 
k}.) It is not hard to see that the fundamental domain is a cube. The action of the 
covering transformations shows you that each face of the cube is glued to the opposite face 
with a one-quarter turn. Eight copies of the cube tile the 3-sphcrc like the eight faces of a 
hypercube. 

The following table shows the results of applying the algorithm to each finite group of 
symmetries of the 2-sphere: 



symmetry 


group of 2-sphere 


spherical 3-manifold 


cyclic 


Z/n 


lens space L(2n, 1) 
fundamental domain is lens-shaped solid 
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dihedral 


Dn 




prism manifold 
fundamental domain is prism with 27z-gon base 


tetrahedral 


(order 


12) 


octahedral space 
fundamental domain is octahedron 
24 of which tile the 3-sphere 
in the pattern of a regular 24-cell 


octahedral 


(order 


24) 


snub cube space 
fundamental domain is snub cube 
48 of which tile the 3-sphere 


dodecahedra] 


(order 


60) 


Poincare dodecahedral space 
fundamental domain is dodecahedron 
120 of which tile the 3-sphere 
in the pattern of a regular 120-cell 



Second category: 

Second category manifolds do not arise from simple left multiplication by groups of 
quaternions. The simplest second category manifolds are generic lens spaces. The group 
of covering transformations of a lens space is cyclic, generated by a single screw motion. 
For the lens space L{p,q), with p and q relatively prime integers satisfying p > q > 0, 
the screw motion translates a distance 27r/p along a given geodesic, while rotating through 
an angle q{2Tr/p) about that same geodesic. The fundamental domain is a lens-shaped 
solid, p copies of which tile the 3-sphere in much the same way that p sectors tile the 
surface of an orange. In the special case that q = 1, the screw motion corresponds to 
left multiplication by the quaternion cos(27r/p) -|- isin(27r/p), so L{p, 1) is a first category 
manifold. Otherwise L{p,q) is a second category manifold. (Warning: L{p,p— 1) is the 
mirror image of L{p, 1) and therefore has the properties of a first category manifold. It 
corresponds to right multiplication by cos(27r/p) sin(27r/p) instead of left multiplication.) 

The remaining second category manifolds are defined by simultaneous left and right 
multiplication by quaternions. That is, the holonomy group of each such manifold is a 
subgroup of a product Hi x H2 x =bl, where Hi and H2 are finite groups of unit length 
quaternions. An element (/ii, ±1) acts according to the rule q — > ^.hiqK^^ . For details. 
Section 4.4 of |22| and Sections 3 and 4 of citenewsphere. 

Flexibility: Spherical 3-manifolds are rigid; they have no flexibility. (As in 2-dimensions 
we make the assumption that spherical geometry refers to a sphere of radius one. Otherwise 
all spherical manifolds would have the flexibility of rescaling - they could be made larger 
or smaller.) 

Symmetry: First category manifolds are always globally homogeneous, that is, each has 
a continuous 3-parameter family of symmetries taken any point to any other point. 

Proposition N.3.3.4 (global homogeneity of first category manifolds). For any two points 
P and Q in a first category manifold M, there is a symmetry of M taking P to Q. 

Proof. Think of M as the quotient of the 3-sphere under the action of a finite group of 
unit length quaternions {gi = 1, g'2) •••) 5n}- The point P has n images represented by the 
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quaternions {gox, ...,gnx}, and the point Q has n images represented by the quaternions 
{goy, gnU}, for some x and y. Right multiphcation by {x~^y) is a symmetry of the 3- 
sphere taking the images of P to the images of Q. More generahy, right multiphcation takes 
every set of equivalent quaternions to some other set of equivalent quaternions, so it defines 
not only a symmetry of the 3-sphere but also a symmetry of the original manifold M. Q.E.D. 

Second category manifolds have smaller continuous families of symmetries. For exam- 
ple, the lens space L(p, q) has a 2-parameter family of continuous symmetries defined by 
translations along, and rotations about, its central geodesic. Both first and second category 
manifolds may have discrete symmetries, which depend on the manifold. 

Dirichlet domain: The Dirichlet domain of a first category manifold never depends 
on the choice of basepoint. Proof: By Proposition N.3.3.4 the manifold has a symmetry 
taking any point to any other point, so the manifold looks the same at all basepoints, and 
the Dirichlet domain construction must yield the same result. By contrast, the Dirichlet 
domain for a second category manifold does depend on the choice of basepoint. 

Observational status: Recent evidence of an approximately flat observable universe 
implies that if the universe is spherical its curvature radius may be comparable to or larger 
than our horizon radius. So even though all spherical manifolds are viable candidates for the 
topology of the universe, the simplest ones, which have large fundamental domains, would be 
much larger than the horizon radius and couldn't be observed directly. A more complicated 
topology, with a larger group of covering transformations and a smaller fundamental domain, 
would be more amenable to direct observation. 



3.2.4 Hyperbolic three-dimensional manifolds 

Manifolds: In 1924 Alexander Friedmann published a paper [|3^ complementing his earlier 
model of an expanding spherical universe with a model of an expanding hyperbolic universe. 
It is a tribute Friedmann's foresight that he explicitly allowed for the possibility of a closed 
hyperbolic universe at a time when no closed hyperbolic 3-manifolds were known. Fortu- 



nately examples were not long in coming. Lobell published the first in 1929 |34] and Seifert 



and Weber published a much simpler example soon thereafter [35|. Nevertheless, during 
the first half of the twentieth century few closed hyperbolic 3-manifolds were know. The 
situation changed dramatically in the mid 1970's with the work of Thurston, who showed 



that most closed 3-manifolds admit a hyperbolic geometry. For an overview, see Ref. [31 
To study known low- volume closed hyperbolic 3-manifolds, see Ref. [ p^ . 

Flexibility: Unlike closed hyperbolic 2-manifolds, closed hyperbolic 3-manifolds are rigid; 
they have no flexibility. As usual we assume curvature radius one, to suppress the trivial 
flexibility of rescaling. 

Symmetry: Closed hyperbolic 3-manifolds never have continuous families of symmetries. 
Truly generic closed hyperbolic 3-manifolds appear to have no discrete symmetries either, 
although the symmetry groups of low- volume manifolds are typically small cyclic or dihedral 
groups. Nevertheless, Kojima has shown that every finite group occurs as the symmetry 
group of a closed hyperbolic 3- manifold |3(;]. 



Dirichlet domain: The Dirichlet domain of a closed hyperbolic 3-manifold always de- 
pends on the choice of basepoint. For carefully chosen basepoints, the Dirichlet domain may 
display some of the manifold's symmetries. The Dirichlet domain may also display sym- 
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metries beyond those of the manifold itself; these so-called hidden symmetries are actually 
symmetries of finite-sheeted covering spaces. 
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Figure 27: Two different surfaces of last scattering, each centered on a different galaxy. 



4 Standard Cosmology and The Cosmic Microwave Back- 
ground 

In standard big bang theory, the universe is created in a hot energetic state, a giant primor- 
dial soup. The plasma is opaque to light as photons scatter off hot charged particles. As the 
universe expands, the soup cools and some 300,000 years into our history, the cosmos cools 
enough so that light no longer scatters efficiently and the universe becomes transparent 
to radiation. This transition, known as decoupling, marks the time when the primordial 
radiation can free stream throughout the universe unimpeded. Decoupling happens quickly 
but not instantaneously and the small spread in time will introduce some complication in 
our attempts to measure topology as discussed in later sections. This ancient radiation 
filling all of space has cooled today to a mere 2.73° K, in the microwave range from which 
it has acquired the name the cosmic microwave background (CMB). As the CMB carries 
information almost from the beginning of time, it is a relic of our deep past. 

The light received at the Earth has traveled the same distance in all directions since last 
scattering and therefore defines a sphere. This sphere is known as the surface of last scatter 
(SLS) and is only defined relative to a given observer. An alien civilization in a neighboring 
cluster of galaxies, if coincidentally also conducting observations of the CMB today, would 



receive photons from its own SLS as shown in fig. 27. This will be particularly relevant to 
the circles method of detecting topology described in section § 6.2.1 . 

The famed CMB carries much coveted information about the universe. It has been 
hoped that all local geometric quantities such as the local curvature, the expansion rate, 
the nature of the dark matter, and of the luminous matter, can be extracted from this one 
source. Clearly though the data will leave the full set of parameters underdetermined and a 
great deal of model dependence is unavoidable. What exactly the CMB can tell us is being 
put to the test with the recent high resolution experiments such as TOCO [^], Boomerang 
and MAXIMA [^8| . These experiments are consistent with an approximately flat observable 
universe. If the universe is hyperbolic, its curvature radius may be comparable to or larger 
than the horizon radius. So even though all hyperbolic manifolds are viable candidates 
for the topology of the universe, they may be too large to see, therein is the challenge. 
The CMB can also be used to decipher the global topology as well as fitting these many 
local parameters. Though statistical measures of topology are inherently model dependent, 
geometric methods are model independent. The pattern based searches of circles in the sky 
and topological pattern formation require no assumptions about the primordial spectrum 
nor the nature of the dark matter. However, these searches will likely be plagued by 
difficulties of their own and it is unclear if these methods will succeed without a great deal 
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of refitting. 

The standard big bang theory also has a powerful successor, the standard inflationary 
theory. The standard inflationary theory suggests there was a period of accelerated ex- 
pansion very early in the universe's history which dilutes all fluctuations, all matter, and 
all energy so that the cosmos suddenly goes from a tumultuous beginning to a nearly flat, 
empty, and biggish universe. Inflation is exited when it becomes energetically favorable for 
the energy driving the expansion to be released into a renewed soup of high energy particles 
and light. This heating of the universe mimics the original primordial soup and the standard 
history of the universe resumes. There are proponents of inflation who would revise the big 
bang still further, essentially removing the moment of genesis altogether so that inflation 
is eternal |39, 4C]. The universe can be viewed as a ginger root of inflationary patches each 
with its own details. As inflation ends in any given patch, the energy is released to heat that 
region. As far as any given observer is concerned, the heating epoch on exit from inflation 
looks like a consequence of a big bang event. Since this is all we can ever observe with the 
CMB our discussion will apply to any of these early universe scripts. 

While theorists still find it tricky to identify a specific working model of inflation and 
to devise a graceful exit to inflation, it is the favored story to tell and some of the reasons 
for this are good. Inflation elegantly resolves the mysterious homogeneity and isotropy of 
our observable universe. Without incredibly precise tuning of the initial conditions, the 
universe would tend to be extremely lumpy and would quickly cool if underdense or would 
quickly collapse if overdense. The observational fact that the universe is billions of years 
old and still just marginally curved, if not outright flat, seems special and so unsettling. 
A more rigorous and compelling argument is made by Guth in his original paper |41] and 
many subsequent articles |^2[. Our homogeneous and isotropic space becomes remarkably 
unlikely, and yet here we are. 

On the other hand lumpiness seems natural since the many regions of the universe 
would not be in causal contact and therefore should have very different local properties 
such as temperature, density etc.. Yet when we look at the CMB we see that in fact the 
temperature appears identical to better than 1 part in 10^ even for regions which would be 
separated far beyond causally connected distances. This would be like finding two ancient 
civilizations on opposite sides of the Earth with nearly identical languages. The civilizations 
must have been in causal contact. Likewise two regions of the universe which seem to have 
equilibrated to precisely the same temperature must have been in causal contact. Inflation 
resolves this situation by taking a small causally connected region and stretching it so large 
that it exceeds the extent of the observable universe today. The CMB is essentially the 
same temperature throughout our patch of the cosmos precisely because it was in causal 
contact. The stretching also naturally renders the universe flat. 

Another motivation for inflation is that, while it makes the universe on average homo- 
geneous and isotropic, it also naturally generates fluctuations about this average. These 
minute perturbations are critical for initiating the collapse of matter into galaxies and clus- 
ters. These seemingly insignificant fiuctuations become the catalysts for all the order and 
structure we observe in the universe. Because of the enormous stretching of scales involved, 
it manages to place fluctuations even on the largest possible scales and therefore gives a 
causal explanation for the seeds of structure formation. (Another mechanism for generating 
the all important initial density fluctuations is via topological defect models although these 
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have lately fallen into disrepute p3| .) 

To emphasize, the beauty is that inflation both explains the average homogeneity and 
isotropy as well as giving a prediction for the deviations from homogeneity, observed both in 
the CMB and in the form of galaxies. It seems a good question to ask , Is topology consistent 
with inflation? Precisely as inflation drives the curvature scale beyond observational reach, 
so too will it drive any topological scale beyond observational reach. The universe may still 
be finite, only we'll never know it. Certainly a model of inflation could be concocted which 
exited just at the critical moment so that the topology of the universe, and presumably 
the curvature too, is just within observation now, when we also happen to be here to 
look. However, such a contraption would be a perversion against the spirit of inflation as 
a restorer of naturalness. If an observable topology scale and inflation are to be conjoined, 
the reasons must be profound and topology must be an integral part of the inflationary 
mechanism. This is conceivable. For instance the cosmic Casimir effect is the contribution 
to the vacuum energy due to any topological boundary. If this effect were to dilute as the 
topology scale expanded then one can envision inflation driven by this contribution to the 
vacuum and ending precisely as the topology scale became some natural big size and the 
vacuum energy became some correspondingly natural small size. To really wrap it all into 
one clean package one could even try to explain the seemingly unnatural Oa = 0.7 this 



way ]lq] (see also Refs. [44, Attempts to use compact topology to aid $7 < 1 inflation 



have also been suggested [^6[. A related mechanism has been put to use in a model with 
topologically compact extra dimensions [47|. Another recent argument to justify flatness 



and local isotropy with topology and without inflation can be found in Ref. |48|. Barrow 
and Kodama found that compact topology can severely restrict the anisotropies which an 
infinite universe will allow. Roughly speaking, the anisotropic modes will not fit in the 
compact space. 

In any case, the inflationary predictions for the statistical distribution of fluctuations 
are often used in fixing the initial fluctuations and so will be relevant here. This assumption 
is a weakness of all statistical constraints on topology and we will discuss these issues in 
due course. The topic for this section is really, regardless of their origin, how do fluctua- 
tions in spacetime translate into temperature fluctuations in the CMB. And once this has 
been established, how does topology alter the standard predictions leaving an archeological 
imprint for us to uncover billions of years later. 

4.1 Standard Cosmological Equations 

This section provides a quick pedestrian review of the CMB and theories for the generation 
of its perturbations. There is nothing specific to topology in this section unless explicitly 
stated. Many detailed reviews on the CMB have been written and the reader is referred to 
these [49, 51, Observing topology in the CMB is discussed in sections 4.4, ^, |6|. 



The gravitational field in general relativity is determined by the local Einstein equations 

G^, = 87rGr^,. (4.1) 

The tensor G^y describes the curvature and evolution of the metric g^^ while T^i, is the 
energy-momentum tensor and accounts for the matter fields. (For thorough discussions on 



general relativity see [m, 54, 53].) It is important to note here that the Einstein equations are 
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local and therefore only determine the curvature of spacetime and do not fix the topology. 
Furthermore the symmetries of g^i, are also local and are nearly always broken in the 
manifold by topological identifications. Today the universe appears to be homogeneous and 
isotropic on the largest-scales. The CMB gives tremendous confirmation of these symmetries 
since the temperature appears identical in every direction to better than 1 part in 10^. 
Homogeneity and isotropy imply that the Earth is not in a privileged position in the cosmos. 
The symmetries of homogeneity and isotropy severely restricts the class of possible solutions 
to the Friedman-Robertson- Walker (FRW) models defined by the gravitational metric 



ds = g^ydx^dx'^ 

dr"^ 



-dr +a\'q) 



1 — Kr^ 



+ r^{d9^ + sin 



(4.2) 



where n = 0, —1, 1 corresponds to a flat, negatively curved (hyperbolic) and positively 
curved (elliptical) space respectively (see fig. |^). Negatively curved space is often termed 
"open" and positively curved space "closed" in reference to their simply connected geome- 
tries. Since we are interested in constructing compact, multiconnected spaces, we will avoid 
the "open", "closed" terminology. 

The overall scale factor, a{rj), describing the expansion of the universe is determined by 
the energy of matter through the remaining Einstein equations. We can always operate in 
comoving units where the manifold is treated as a static constant curvature manifold and 
all of the dynamics is hidden in the conformal scale factor a(r/). An alternative expression 
for the metric is 

ds'^ = a^i-q) [-d7?2 + da^] (4.3) 
with the spatial part of the metric 

da^ = dx^ + f{x){d9^ + sin^ ^#^) (4-4) 
where we have used conformal time d-q = dt/a{t) and 

{X^ r = X flat 

sinh^ X r = sinhx hyperbolic (4.5) 
sin^ X r = sin x spherical 

See also Appendix |^. 

The dynamical evolution of the space is given by the Einstein equation determining the 
scale factor, 

+ K/a^ = ^p. (4.6) 

The different curvatures correspond to different values of the global energy density. Tradi- 
tionally, Vt = p/ pc '\s defined with pc the critical density required to render the universe flat, 
so that > 1 corresponds to k = 1, J7 = 1 corresponds to k = 0, and < 1 corresponds to 
K = — 1. The curvature radius is Rcnw = a/\n\^^'^ ■ From eqn (|4.6| ), SirGpc/'i = so that 
H^a'^{Q — 1) = K and the curvature radius can be be expressed in terms of Q as 

-^curv ! , 1 /n • (4' 7) 
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When working in comoving units we take the comoving curvature radius to be unity for 
curved space or oo for flat space. 
Conservation of energy requires 

p + ?,H{p + p) = {) (4.8) 

for the perfect fluid energy- momentum tensor = {p,p,p,p)- Decoupling occurs during 
the matter dominated error for which p = and so p cc by eqn (^^) . The solution for 
a (77) during matter domination is 




0(77) oc < r/72 K = (4.9) 



with the conformal Hubble expansion 7i = a' /a = aH, 

{sinhT//(coshr/ — 1) k = —1 
2/r] K = . (4.10) 

sinr//(l — cos 77) k = 1 

In comoving units, the universe is static. In physical units, the universe expands and 
the cosmic background radiation redshifts. The physical wavelength is shifted according to 

m ait) 



^initial o(iinitial) 

from which the cosmological redshift is defined: 

l + z = ^. (4.12) 

a{t) 

A subscript will always be used to denote values today. The value of r] during matter 
domination can be written in terms of and the redshift using r] = J dt/a = J da/ (da) 



and eqn (4.6): 



arccosh (1 + g^) no<l 



(1+2)1/2 



^0 = 1 (4.13) 



cos-ifl-^fei)) 17n>l 



no(i+z)) "0 



If curvature or topology are to be observable, the geometric scales must be smaller than 
the diameter of the SLS. In flat space, there is no natural scale since i?curv = 00 and it 
would seem a random coincidence if the topology scale just happened to be observable. In 
curved space, one might expect the curvature and topology scales to be comparable. This 
may not make it any more natural to observe geometric effects, but it would mean that if 
curvature is observable, then topology may be also. 

For the purpose of observational topology it is useful to list some relevant scales. The 
radius of the surface of last scatter is defined as the the distance light travels between the 
time of decoupling and today, 

/■*« dt 

D, = a{t) — (4.14) 
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The distance depends on the curvature and evolution of the universe through eqn (4.9). In 
comoving units the photon travels a distance 




'*o dt 



fVO 

/ dr] = Ai] 



(4.15) 



a{t) 



td «(*) 



with Ar] = r]Q — r](i. The diameter of the SLS, 2Ar], essentially defines the extent of the 
observable universe. A loose criterion that can be used to gauge if topology will influence 
the CMB is that the in-radius be less than Aij. 



The volume enclosed by the SLS is the integral over Vsls = J yj —gdrdOdcj). For a radius 
of Arj, 



The number of clones of the fundamental domain that can be observed in a small compact 
cosmos can be estimated by the number of copies that can fit within the SLS: Vsls/^m- 
We will often refer back to these scale comparisons in the coming discussions. 

4.2 Fluctuations in the CMB 

Regardless of the origin, there are small deviations from homogeneity and isotropy which 
are treated by perturbing the metric and matter tensors about the FRW solutions. The 
equations of motion for these perturbations to linear order are determined by 5Gfj_y = 
SnGSTfj^iy . While the equations are extremely complicated there is one quantity, the gauge 
invariant potential which is a coordinate invariant combination of components of dg^u 
and is of central importance for our considerations. It is so named since in the Newtonian 
limit it corresponds to the usual gravitational potential. 

The hot and cold spots in the CMB are caused by these perturbations in the gravita- 
tional potential. The fluctuations Doppler shift the photons generating an anisotropy both 
at decoupling and as light traverses the time changing gravitational field. There are many 
carefully derived results for the temperature fluctuations in the literature based on rela- 
tivistic kinetic theory and relativistic perturbation theory ||57| , Since these topics 

comprise multiple review papers on their own, we will not go through the detailed deriva- 
tions but will instead try to provide a cohesive, intuitive motivation for each of the relevant 
concepts. 

To derive the famous Sachs- Wolfe effect we followed the pedagogical discussion of 
White and Hu p^]. Photons gain energy and are therefore blueshifted as they fall into a 
potential well and lose energy and are therefore redshifted as they climb out. The Doppler 
shifts will cancel in a static homogeneous and isotropic space. However at the time that the 
photons last scatter, some photons will be in potential wells and some will be in potential 
peaks. Therefore this snapshot of the metric fluctuations becomes frozen into the CMB. It 
is this snapshot or fossil record we observe. Energy conservation gives a Doppler shift of 



Vsls 




(4.16) 



5T 



5T 



(4.17) 




T 
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The local gravitational potential makes an isotropic contribution to 6T/T and will be 
left off hereafter. By adiabaticity, aT =constant, it follows that 



f =-^. (4.18) 

T i a 

Assuming flat space for simplicity, a oc t"^^'^ during matter domination. The shift in a is 
then 

5a 2 6t , , 

- = i;-- 4.19 

a 6 t 

The shift 6t can be understood as a shift relative to the cosmic time due to the perturbed 

gravitational potential. In a gravitational potential clocks appear to run slow by 

y ~ $ (4.20) 



from which it follows that 

^. (4.21) 



6T 



This term is the surface Sachs- Wolfe effect from the time of decoupling. There is an ad- 
ditional contribution to the temperature fluctuation known as the integrated Sachs- Wolfe 
(ISW) effect which is accumulated as the photon traverses space through a decaying poten- 
tial. Including this latter contribution the full Sachs- Wolfe effect is 



5T{n) 1 



$(r?on) +2 r° (ir/$'(7?,x) (4.22) 



where a prime denotes differentiation with respect to conformal time ||58[. In flat space 
with ordinary energy, the only contribution to 5T /T is from the surface Sachs- Wolfe effect 
(the first term in eqn. [4.221) . The ISW (the second term in eqn. 4.22 ) makes a significant 



contribution to large-scale fluctuations if space is curved and/or when a cosmological con- 
stant dominates the energy density. In the absence of a significant cosmological constant, 
the ISW is important as perturbations decay along the line of sight during the curvature 
dominated epoch at 

l + z~(l-f]o)M (4.23) 

for underdense cosmologies. This will prove to be important to the issue of topology since 
the ISW may camouflage the conspicuous marks of topology imprinted in the surface Sachs- 
Wolfe effect. 

Therefore the temperature fluctuations observed today can be predicted from the pri- 
mordial fluctuations in the metric. The shape of <I>(?7,x) will depend on the geometry of 
space as well as some initial spectrum. The perturbed Einstein equations give the equation 
of motion 

^" + m{i + cl)^' - clv^^ + {27i' + (1 + cDn^)^ = (4.24) 

where 7i = a' /a and Cg is the speed of sound in the cosmological fluid. In a matter dominated 
era Cg = 0. The potential is separable and can be written $ = F{r])'^{x) with F{r]) the 
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solution to eqn (g]2|). The ^'(x) can always be expanded in terms of the eigenmodes to 
the Laplacian; that is, 

$(?7,x) = F(r/) J cPk Vfc(^) (4.25) 
with the eigenmodes if^k satisfying 

VVfc = -feVfc- (4.26) 

The Laplacian = g^^D^D^ depends on the curvature through the covariant derivatives 
D^. For now we assume all the /c's are continuous eigenfunctions although this will not be 
the case when the space is topologically identified. The ^>fc are initial amplitudes given by 
the statistical profile of the initial fluctuation spectrum. All assumptions and/or predictions 
for the initial perturbations are contained in the 

During inflation quantum fluctuations in the field are amplified by the accelerated ex- 
pansion. These fluctuations about the ground state of the field theory are known to be 
Gaussian distributed. The amplitude of the fluctuations are related to the speciflc infla- 
tionary model but for the most part are independent of the scale k. So inflation predicts a 
Gaussian distribution of fluctuations independent of scale. Speciflcally this means the $g 
are drawn from a random Gaussian ensemble consistent with 

^^^)^V^{k)5\k-k') . (4.27) 

The angular bracket <> denote an ensemble average and is the predicted power spec- 
trum. De Sitter inflation delivers a flat, Harrison-Zeldovich spectrum which corresponds 
to =constant. Other forms for oc have been derived from speciflc inflationary 
models with the spectral index n ^ 1 and give tilted as opposed to flat spectra. 

From the preceding we can fully predict $ and therefore 6T/T given the curvature of 
spacetime and the initial spectrum. Armed with this prediction, fluctuations on a given 
manifold can be compared to the data. For data comparison it is customary to consider 
the correlation function. Although 5T{x.) /T permeates space, we only observe fluctuations 
from our SLS at location x = Ar/n with n a unit directional vector. Since the fluctuations 
are taken to be Gaussian, the correlation function contains all of the information about 
the temperature fluctuations. The ensemble average correlation function between any two 
points on the sky is 

C{h,h') = ^^(n)^(n')). (4.28) 

The theoretically predicted C{n,h') can then be statistically compared to the data to esti- 
mate the likelihood a given model is responsible for the world we live in. 

Because fluctuations are observed on a sphere it is customary to decompose the data 
into spherical harmonics, 

ST 

— (n) = ^ atmYim{n) (4.29) 

where the l^m's are the usual spherical harmonics. The spherical harmonics form a complete 
set of states on the sphere and are orthogonal so that 

sin9deYl,^,{d,<P)Y,m{e,(t)) = 5ei5^>m (4.30) 
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and unit normalized. Using the orthogonality of the l^mj invert (|4.2g| ) to find 

/ST 

dnYe_m{n) — ih) . (4.31) 

If the probability distribution is Gaussian, then the Fourier correlation function 

Ce = (M^) (4.32) 



contains complete information about the fluctuations. The <> again denotes an ensemble 
average. Either non-Gaussianity or the breaking of homogeneity and isotropy with topo- 
logical identifications will mean that the C^'s do not provide complete information. We will 
emphasize this in section §|^. 

Assuming for now a simply connected topology, homogeneity and isotropy of the metric 
allows one to perform an angular average which is in effect an average over m's without loss 
of information. Such an average gives an estimator for the ensemble average angular power 
spectrum Ci, 

m=£ 

Ci= |a£mlV(2Z + l). (4.33) 

m=—£ 



Using equation ( 4.31 ) gives 



Q = 7^^E[/ dn'Ye-mih') J dOy,^(n)^^(n)^(n')) 



(4.34) 



The parameter i can be associated with the angular size of a given fluctuation, with ^ ~ vr/^. 
Large angle fluctuations are associated with low i and small angle fluctuations are associated 
with high i. The lower the value of £, the fewer contributions there are to the sum ( 4.33| ). 



As a result, the estimator of the true ensemble average is poorer for low £ than it is for high 
i where there are many contributions to the average. This is known as cosmic variance and 
can be included as an error bar. For a homogeneous, isotropic space cosmic variance can 
be estimated as 



Q^2/(2£ + l). (4.35) 

We only have one universe available to measure and we cannot be sure that our universe is 
not just a randomly large deviation from average. It is difficult to interpret the significance 
of low i observations such as the lack of power in the quadrupole as observed by COBE— dmr . 
The ambiguity of cosmic variance is worsened with topological identifications as emphasized 



in HI. 



These are all the tools needed to determine 6T/T for a specific theoretical model, the 
correlation function C(n, h'), and the angular average power spectrum Ci. These predictions 
can then be compared with data. When we come to the influence of topology we will see 
the increasing importance of using full sky maps of 5T/T and the full correlation function 
C{n,h'). 

We have described the scalar modes above. There are also tensor modes which are 
gravitational waves and vector modes which are rotational perturbations and do not grow 
with time. The scalar modes are the only ones which couple to the energy density and 
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pressure and we restrict ourselves to scalar modes hereafter. Topology may in fact influence 
the gravitational wave background but to date this remains unexplored territory. 

The primordial fluctuations catalyze the formation of galaxies and galaxy clusters. So 
the minute quantum fluctuations are amplified into the gigantic structures we see today. 
The primordial spectrum can be tested by measuring the fluctuations in the CMB but there 
are other manifestations. For instance, the eventual nonlinear growth of perturbations can 
be simulated numerically to test theoretical predictions for structure formation against 
astronomical observations. It would be interesting to know if topological features could 
sculpt the distribution of structure as well |]5^ . 

In the following subsections we find the eigenmode decomposition in simply connected 
spaces. Perpetuating a cruel prejudice against S'^ we consider only and E'^. However, 
the authors of Ref. catalog all of the topologies and discuss detection strategies based 
on the crystallographic methods. The motivation for neglecting topologies in connection 
with the CMB is observational as there is no outstanding evidence which supports an > 1 
universe |61]. The current debate is whether Q = 1 or whether Q < 1. The future satellites 



which aim to refine measurements of the CMB intend to resolve this debate. The overriding 
theoretical prejudice is for = 1 consistent with infiation. More recently the supernova 



data has drummed up enthusiasm for I^a = 0.7 and = 0.3 so that the total Q = 1 |64|. 
The recent small angle experiments Boomerang and MAXIMA corroborate these values if 
certain prior assumptions constrain the statistical analysis. More quantitatively, the current 
numbers are O = 0.90 ± 0.15, ^bhl = 0.025 ± 0.010 ^cDMhl = 0.13 ± 0.10 and a spectral 



index of n = 0.99 it 0.09 [65| where the subscript b denotes the baryonic contribution and 
the subscript CDM denotes the net contribution from cold dark matter. Still, the data are 
clearly underdetermined and it is an unresolved issue how many different models can match 
the data. 

In Ref. [^] an argument is made for a constraint on JIq with some model-independence. 
They argue that a constraint is imposed if one requires that a local maximum detected 
in the correlation function of large scale structure (~ 100 — 200/i~^Mpc) occurs at the 
same comoving positions at different redshifts. The argument is independent of both CMB 
and supernova la data, and favors a hyperbolic universe. Their result is a cosmology with 
Jlo = 0.9 lb 0.15 (95% confidence) (with Omegap^ in the vicinity of ~ 0.65). 

Instead of entering further into the parameter estimation debate, we will review the 
recent investigations in cosmic topology which have focused on the fiat and the hyperbolic 
manifolds. For an interesting catalog of topologies for more general Bianchi classes see Ref. 
11 



4.3 Observing the CMB 

Many billions of years after last scattering, we build the COsmic Background Explorer 
(COBE) satellite to confirm the earlier detections by Penzias and Wilson of the microwave 
background radiation. The all sky map generated by the COBE satellite confirmed that the 
average temperature was homogeneous and isotropic at 2.728 ° K and that the spectrum was 
extremely thermal, a result that drew spontaneous applause when presented in 1992. The 
data also confirmed the theoretical expectation that there are in fact minute fluctuations as 
predicted by infiation at the 10~^ level that appear to be scale invariant. COBE measures 
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fluctuations on large scales, £ < 30. Low is correspond to fluctuations on scales far outside 
the horizon at the time of decoupling and therefore can only be due to fluctuations in the 
metric and density. So COBE which measures ^ < 30 is a probe of the largest geometric 
features and therefore is quite important as a probe of topology as well. Causal microphysics 
becomes important for i > 100. At these large £ our calculation of 6T/T in eqn. (|4.22| ) is 
insufficient as it neglects microphysical effects and a more detailed analysis is required. Of 
particular importance are the high i Doppler peaks. Before recombination, sound waves in 
the baryon-photon fluid induce additional Doppler shifts which produce a peak in power on 
a scale related to the size of the horizon at the time of decoupling. The height of the peaks 
depends on the speciflc model parameters. The location of the flrst peak depends primarily 

— 1/2 

on curvature, ipeak ~ 

2201^0 ' , with some small Hq dependence. COBE— DMRdoes not 
measure high £ fluctuations but two important all-sky satellites will. Microwave Anisotropy 
probe (MAP) and Planck Surveyor. The recent balloon borne experiments locate £ ~ 200 
and therefore contribute more evidence in favor of f^o ~ 1- 

The task of the Differential Microwave Radiometer (DMR) experiment on the COBE 
satellite was to measure the large angle anisotropy of the entire sky ^ij |68|| . The tempera- 
ture COBE-DMR measures in a given pixel can be written 

ST \ 

-Tfr) =Y1 ^imBiYimixi) + Hi (4.36) 

where is the experimental beam pattern and the noise in each pixel is n^. The DMR 
horns are characterized by an imperfect Gaussian beam pattern. The noise is assumed to 
be uncorrelated and Gaussian with mean < rii >= and variance < riirij >= a5ij. 

The experiment consists of three pairs of antennae. Each pair measures the temperature 
difference in two directions separated by 60°. Each antenna has a 7° beam and the data is 
smoothed on 10°. A full sky scanned was performed several times over 4 years. The COBE— 



DMR data is provided in the form of 6 maps in 3 frequencies, 31, 53 and 90 GHz |68, 67 1. 
Each map has Np = 6144 of size (2.6")^. Compressing to resolution 5 pixels, Np = 1536 
pixels of size (5.2°)^, loses no information and is sometimes implemented in data analysis. 
The observed C^s have been determined from the COBE-DMRdata by Gorski |p9|, by 



Tegmark |7C], and by Bond, Jaffe and Knox [^]. The monopole {£ = 0) is just the average 
temperature itself and can be discarded from the data. The largest scale anisotropy is the 
dipole {£ = 1) generated by our solar system's peculiar velocity relative to the CMB. Since 
this is not cosmological in origin, the entire dipole is discarded from the data. The flrst 
relevant cosmological observations begin with i = 2, the quadrupole, and it is curious to 
note that the measured quadrupole is low. This may just be cosmic variance in action but 
compact manifolds do happen to predict low power on large scales. This is discussed at 
length in section §|5[ The six maps can be compressed into one weighted-sum map with 
the monopole and dipole subtracted. A galactic cut of a region ±20° around the Galactic 
plane is also needed to eliminate Galactic emission. This finally messaged data is ready for 
comparison to theoretical predictions. 

The signal-to-noise of COBE-dmrIs ^ = 2. By comparison, the projected sensitivity of 
MAP is = 15 with a resolution of 0.5°, a tenfold improvement in signal-to-noise and 30 
in resolution. Planck has even higher resolution but will launch much later. 
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Figure 28: The three panels are, from top to bottom, the combined data before the 
dipole is subtracted, the data after dipole subtraction but before the galactic cut, the 
map minus both dipole and galactic emission. Taken from the COBE webpage at 
http: / /space.gsfc.nasa.gov/astro/cobe/ 
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4.4 Topology and the CMB 

Several aspects of the fluctuation spectrum are altered when the manifold is multiconnected. 
The most conspicuous and consistent signatures of topological identifications are as follows: 
(1) Multiconnectedness destroys global isotropy for all but the projective space and de- 
stroys global homogeneity for all but the the projective space and the hypertorus. (2) The 
spectrum of fluctuations is discrete reflecting the natural harmonics of the flnite space. (3) 
Since the fluctuations must fit within the finite space, there is a cutoff in perturbations with 
wavelengths which exceeds the topological scale in a given direction. (4) Geometric patterns 
are encrypted in the spatial correlations. The patterns reflect the repeated occurrence of 
topologically lensed hot and cold spots. 

The proposed methods of scanning the CMB for evidence of topology can be split 
into direct statistical methods and geometric methods. The direct methods begin with a 
theoretical model, compute the temperature fluctuations as was done in §5.1.1 and § |6.1.1| for 



the simply connected spaces, and determine the likelihood of the model parameters against 
the data. There are shortcomings with this brute force approach. All conclusions are model 
dependent and there are an infinite number of models. The model is not just the manifold 
but also the orientation and the location of the Earth and the initial primordial spectrum. 
While some conclusions can still be drawn there is an appeal to a model independent attempt 
at observing topology. That is where the geometric methods come in. The geometric 
methods treat topological lensing much like gravitational lensing. One just observes the 
lensed images and their distribution to reconstruct the geometry of the intervening lens. 



The most amusing example is the correlated circles of §6.2.1 although other patterns can 



emerge as well as described in §5.2.2 



Before preceding to review the known approaches to observing topology it is worth 
expanding on point (1) above. The global anisotropy and inhomogeneity means that the 
correlation function depends on the location of the observer and the orientation of the 
manifold. Additionally, C{n,fi') depends on both n and fi' and not just the angle between 
them. As a result, the angular average performed in the definition of the multipole moments 

discards important topological information. The information lost in the angular average 



can be quantified by an enhanced cosmic variance (C|) [18|. The correlation function can 
be decomposed into isotropic and anisotropic pieces. 



By isotropy 



and by anisotropy 



C{h, h') = C\h, n') + C^(n, n'). (4.37) 



C\n, h') = Y: l^^fiPdn • n') (4.38) 



i 



C^{h,n) = j dO-n j dnn'C^{h,n)Pi{n ■ n) = (4.39) 

so that the anisotropic piece is orthogonal to the Legendre polynomials. The expectation 
value of the estimator for the Cg depends solely on the isotropic piece by construction 

i^e) = ^^^^ / / dnn,C{h, h')P,{h ■ h'), (4.40) 
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but the variance 

var(Q) = Icj) - (cA' (4.41) 



contains anisotropic pieces. This can be interpreted as very large error bars due to cosmic 



variance. As a resuh, conclusions based on the C^'s alone are weak [17, 18|. Many of the 
statistical analyses described below do only examine the C^s and so can only rule out a 
topology but never really confirm topology. 

However, others have argued that this increased error is not very large for compact 
hyperbolic spaces. An argument by Inoue for instance expresses the variance as the sum 
of a geometric variance and the usual cosmic variance. The geometric variance contains 
the additional variance due to topology and is intrinsically small [|6^]. The smallness of 
the geometric variance can be traced to the randomness of the eigenmodes on compact 
hyperbolic models (explored further in § |6.1.3| ) which generates what Inoue refers to as 
geometric Gaussianity ||6^. 
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5 Observing flat topologies in the CMB 



Since the flat spaces are most easily constrained using the COBE-DMRdata we begin with 
them and discuss hyperbolic manifolds in all their glory separately. While small universes 
can be ruled out, it is rather fascinating to note that large cases are marginally consistent 
with the data. After all, the observed quadrupole is low and so the infrared truncation in 
power seen in flat spaces and described below can actually create a better fit to the data. 
However, such a marginal flat space, just coinciding with the size of the observable universe 
is unaesthetic at best. In fact, there is no natural scale for the size of a flat universe and 
they are not the favored small universe candidates for this reason. Still, they provide an 
important and accessible testing ground for methods of observation and we discuss them 
next. 



5.1 Direct methods in flat space 

5.1.1 Simply connected fiat space 

In flat, simply connected space, the spectrum of fluctuations is well known. We need to find 
the elements of eqn. (|4.25 ) with F{r]) the solution to eqn. (|42^ ). In flat space 7i = 2/r] and 



the solution to eqn ( [4.24 ) decays as ^' oc r]~^. Consequently, in flat space F(t]) is effectively 
constant and there is no ISW effect during radiation or matter domination. (There is 
however an ISW effect if the universe is dominated by a cosmological constant. This may 
turn out to be important in salvaging finite flat models if the universe is accelerating as the 
recent supernovae observations indicate The Laplacian in flat space is simply 



[dl + d^y + a^) = -feVfc (5.42) 

with solutions 



ipj: = exp [ik ■ xj . (5.43) 
The potential of eqn. (4.25) can then be expanded in terms of these as 



$(x) = /"°° d^/cl>^exp (iA?-x) (5.44) 

All of the assumptions about the statistics and shape of the spectrum are contained in 
the <I>g. All other quantities are determined by the geometry of the space. In accordance 
with the inflationary prediction the cE>g are assigned an independent Gaussian probability 
distribution consistent with the flat space normalization 

hn) = ^'P^(^)^'ik-k') (5.45) 



11 



Prom eqn. (4.22) the Sachs- Wolfe effect is simply 

ST{h) $(n) 



(5.46) 
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The correlation function between any two points on the SLS is then, from eqn ( fl.28| ), 

-p-'Pcj) exp yiArjk ■ {h — n)j (5-47) 

up to a normahzation. For a homogeneous and isotropic space, the correlation function 
depends only on the angular separation between n and h' . Consequently all the information 
on the theoretical sky is in the angular average. Using the orthogonality relations of the 
Legendre polynomials in eqn ( [4.34| ) gives 



If we assume the initial fluctuation is a powerlaw oc A;" , this can be integrated to give 

r(£+(n-l)/2)r((9-n)/2) 
°^ r(^ + (5 - n)/2) r((3 + n)/2) ^^'^^^ 



As described in §4.2, the COBE— DMRdata has been analyzed over the years by many 
groups to determine the C^'s observed. Statistical likelihood analyses are then performed 
to compare the theoretical model with the COBE-dmr data. The infinite, fiat model is 
generally found to be consistent with the data for a nearly fiat power spectrum (n ~ 1). 
The question addressed in the next section, is how well finite flat models match the data. 
To flash forward, the answer is essentially that finite flat models match the data well if 
they are comparable in size to half the observable universe although some studies have put 
even stronger limits as described below. All limits are subject to caveats based on the 
assumptions made about the perturbation spectrum - none of which have yet been tested 
observationally. 

5.1.2 Compact, fiat spaces 

We have already reviewed the standard decomposition of the temperature fluctuations in a 
simply connected space. When the manifold is compact and the SLS exceeds the dimensions 
of the space, then the global topology is reflected in the CMB sky. There are three notable 
alterations to the predicted fluctuations when the manifold is compact: (1) The eigenvalue 
spectrum is discrete not continuous. (2) There is a cutoff in the power of fluctuations 
on wavelengths which exceed the natural size of the space. (3) The correlation function 
C(n, n') depends on orientation and so depends on both fi and h' explicitly. It is not simply 
a function of the angular separation. 

The exact eigenmodes can be found for all of the 6 compact orientable flat spaces |]73| , 
74 1 . This allows a direct statistical comparison of the theoretical predictions for compact, 



orientable spaces with the COBE-DMRdata. With any statistical comparison of a model 
with the data, the conclusions are model dependent. The reader should bear in mind that 
the bounds quoted assume equal-sided spaces, = 0, and a flat Gaussian distributed 
power spectrum. 
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The earliest bounds on the hypertorus constrained the topology scale to be ^ O.SAr] 
1 75 1 which is still less than the diameter of the SLS. There could still be as many as eight 
copies of our universe within the observable horizon. The analysis was later extended to the 
other compact, orientable flat spaces where similar bounds were obtained 74 1. Stronger 



bounds were put on the hypertorus by comparing the full covariance matrix against the 
2-year COBE-DMRdata ||76|. The length of a side was set to be > 1.2Ari. 



The most conspicuous feature in the fluctuation spectrum is a suppression of power on 
large scales since large fluctuations cannot fit into the finite box. Such an infrared cutoff 



can be deduced from Cheeger's inequality [77| 



T' = ¥ min(y(Mi),y(M,)) ' ^^'^^^ 



with the infimum taken over all possible surfaces S that divide the space into two subspaces 
Ail and Ai2 where M. = Aii Ci M.2- S is the boundary of the two subspaces, S = dMi = 
dM2- The isoperimetric constant he depends on the geometry more than the topology. 
Intuitively speaking, in very long thin manifolds he can be quite small leading to a lowered 
spectrum of eigenvalues. However, this is highly correlated with thin bottleneck structures 



and is not a common feature of spaces with a more regular shape [77, 78|. From Cheeger's 



inequality it follows that all flat hypertori have kmin > 2/L with L the longest side of the 



torus as argued in |17|. 



Although the cutoff would seem a good indicator of topology, it just so happens that 
the longest wavelength fluctuation observed, namely the quadrupole, is in fact low. Some 



might even take this as evidence for topology [79|. Cosmic variance is also large on large 
scales. Consequently, a fundamental domain the size of the observable universe is actually 
consistent with the observed COBE— dmrC^s |7^. However, since the fundamental domain 
has a particular orientation on the sky, the correlation is not simply a function of the angular 
separation between h and h' as it is in the infinite case. Therefore conclusions based on the 
C^s alone are weaker than a statistical analysis based on the full C{h,h'). The strongest 
bounds on the hypertorus were placed by comparing the full correlation function C(n, h') to 



the data. It was found that L > 1.3Ar/ |7q, 17, which is still less than the diameter of 



the observable universe, ~ 2A?7. Asymmetric spaces were constrained by de Oliviera-Costa 



et. al. as described in section S5.2 



Despite the bounds placed on the equal-sided compact, flat spaces, they still provide an 
excellent testing ground for geometric measures of topology. They may also be saved by 
Q\ 7^ models. We compile a list of all the eigenmodes, eigenvalues and relations for the 



compact flat spaces for completeness. These solutions are taken from Ref. |74]. The spectra 
were also found in Ref. |^] for the other compact spaces, however, the critical relations 
were overlooked leading to confusion about the cutoff in the long wavelength modes. It 
was mistakenly concluded because of the missing relations that in the twisted spaces longer 
wavelengths could fit in the fundamental domain since a wave needs to wrap more than 
once before coming back to a fully periodic identification |l^. However it was shown in 
Ref. [74|, that these long modes are forbidden by the relations among the and all of 



the compact topologies have roughly the same long wavelength cutoff. 

As a result of the global topology, all of these spaces are anisotropic and all except for the 
hypertorus are inhomogeneous. Topology can be implemented by imposing the boundary 
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Figure 29: Tiling flat space with parallelepipeds. 



conditions 

$^(x) = $^(5x) v^er. (5.51) 

Compact topology always restricts the eigenvalues to a discrete spectrum: 

^(^) = EV'" • (5-52) 

k 

In additional, relations are imposed on the 

Demonstrating the discretization explicitly for the hypertorus, the identifications of 
the cube are expressed in the three boundary conditions ^{x,y,z) = ^{x + L,j.,y,z) = 
$(x,y + Ly,z) = ^{x,y,z + L^). Imposing the first boundary condition on of eqn. 
(|5.52| ) gives e~^^''^ = ^-^^^ix+h) ^\^if.\^ requires kx = j^^x with Ux an integer. The other 
two boundary conditions provide the discrete spectrum 

27r , 27r , 27r 

ly 



kx — — nx ky — — Uy kz — — Tiz (5.53) 

'.r '?/ I'Z 



with the Ui running over all integers [pO, SI, U^. It is clear that there is a minimum 



eigenvalue and hence a maximum wavelength which can fit inside the fundamental domain 
defined by the parallelepiped [^] : 



fcmin = 27r min (^J-^ ^max = max (Lj 



Although global anisotropy is broken by topology, we can still form the C^s to make a com- 
parison, if incomplete, with the COBE— dmrC^s. Expanding the exponential and Legendre 
polynomials in terms of spherical harmonics, the Ci eqn. ( [4.34 ) becomes 



Ce^Y.Y.'^j'(^vk) (5.54) 

k ¥ 

Three other spacetimes are constructed from a parallelepiped. The first twisted paral- 
lelepiped has opposite faces identified with one pair of faces twisted by vr before gluing. The 
eigenmodes are k = 2iT{j /h,w /h,n/2c), with the additional relation ^j^n = ^-j-wn e.™"^. 
The minimum mode consistent with the relations is /cmin = 27r/L. The C^s become 

<^^°^Z^-p 21 > I 1^ I^Am(k)| (1 + e ^ >). (5.55) 

jwn m=—£ 
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Figure 30: Tiling flat space with hexagonal prisms. 

A parallelepiped with one face rotated by 7r/2 has discrete eigenmodes k = 2Tr{j/h, w/b^n/Ac) 
with the additional relations ^jwn = ^w-jne''"'^^^ = ^-w-jn.e^"''^ = ^-wjne'^"''^^^ ■ Again 
the minimum mode is k^in = 27r/L. The C^s are given by 

k It -\- \ 

The last parallelepiped has a fundamental domain of volume 2hhc and is thus a double 



parallelepiped. The identification rules involve three rotations through vr as shown in |19|. 
They are so that (x, y, z) {x + h, —y, —z), {x, y, z) {—x, y + 6, — (z + c)), and {x, y, z) — > 
(— (x + h), —{y + b), z + c). The resultant discrete spectrum is A; = TT{j /h,w/b,n/4:c) with 
relations ^jwn = ^j-w-n e*'^-' = 4>_j^_„ e^'^('^+") = ^-j-wn e*'^(-^+"'+"'). The minimum 
mode is kmm = V^tt/L. Since the volume of the space is 2L^, the angular averaged long 
wavelength is \/2/2^/'^ — 1.1 times the length of the fundamental domain. The angular 
power spectrum is 

Ct ^ -p ^ ^ y^,„(k) X (5.57) 

jwn m=—£ 

The last two possible compact flat spaces are based on a hexagonal tiling. With the 
opposite sides of the hexagon identifled and the prism faces rotated relative to each other 
by 27r/3, the potential can be written 



^ = 5Z inananze^^'^'exp n2 (^x - -^y^ + ns (^x + -^y 



(5.58) 



with the eigenmodes k = 27r((n2 + n^)/h, (— n2 + n^)/b^ Uz/'ic). The relations on this space 
are I>n2,n3,n, = $„3,_(„2+n3),n.e*2''"'^^ = ^-(n2+n3),n2,nze'^™'^^ and lead to 

Cio, ^iflT E |l^.,n.(k)P(l + e^^(-+")-/3 + e^^(-+-W3). (5.59) 

Lastly, the prism faces are glued after rotation by vr/S. The potential can still be 
written as ( p. 58 ), with eigenmodes k = 27r((n2 + n^)/h, (—712 + n^) /b,nz/Qc) and a set of 
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relations '^n2,nz,n^ = *P(,22+n3),-n2,nz^ ^ ~ '^n3,-(n2+n3),nz6 ^ ~ *P-n2,-n3,"z ^ ^ ~ 

i>_(„2+„3),„2,„,e^4™-/3 ^ l>_„3^(„2+„3),nze*^™^/^. The Qs are given by 

E -p E (5.60) 

n2n3nz m=-i 
^ -j^ _|_ gi(nz+m)7r/3 _|_ gi2(nz+m)7r/3 _|_ ^i(n^+m)-n _|_ gj4{nz+m)7r/3 _|_ gj5(nz+m)7r/3-j 

The volume of both of these topologies is 

In summary, both the vr- twisted and 7r/2- twisted tori have C^s that are almost identical 
to that of the torus. The harmonics do show distinctions between the parallelepiped topolo- 
gies, but unfortunately the variations fall well within cosmic variances. In addition to the 
damping at low harmonics of the discrete spectrum create both dips and enhancements 
extending up to high values of I. The dips are a natural consequence of the discretization 
and the resultant absence of certain modes in the spectrum. The enhancements are due to 
the distribution of multiple images dictated by the geometry of the fundamental domain. 

As a result of the low observed quadrupole, the predicted cutoff alone is not enough 
to rule out compact, flat models. The likelihood of the in the full COBE—dmr range 



was compared to the relative likelihood of a flat, infinite cosmology in Ref. |75, 74 1. 
Compact, flat spaces are tens of times less likely than their infinite counterparts if the 
topology scale exceeds about half the diameter of the SLS. Since compact topologies do not 
give isotropic temperature fluctuations, this lends ambiguity to any likelihood analysis. The 
conclusions drawn were quite conservative, ruling out (. > 0.8Ar/; that is, the length of a 
side must be ^ 0.4 the diameter of the observable universe. Again, it should be emphasized 
that these conclusions are contingent on the correctness of the assumptions about the initial 
perturbation spectrum. 

While the angular power spectrum is sufficient to constrain symmetric, flat topology it 
is in general a poor discriminant. The average over the sky fails to recognize the strong 
inhomogeneity and anisotropy manifest in these cosmologies. Direct attacks on anisotropic 
spaces inspires more geometric approaches as discussed next. 

5.2 Geometric methods in flat space 



One of the first geometric, or pattern driven searches for topology was initiated in Ref. |82|. 
For an anisotropic hypertorus, a symmetry plane or a symmetry axis can be identified in 
the CMB |4^, |8^, For these asymmetric spaces, a search for patterns was emphasized 
in Ref. They develop a statistic that is independent of cell orientations but is sensitive 
to the plane and axis symmetries of the rectangular spaces. For anisotropic models there 
is a stronger dependence on the orientation of the manifold relative to the Galaxy cut. 
The statistic S{hi) searches for reflection symmetries in a plane perpendicular to hi and is 
defined by 

^^^'^ ^ ^ § cT2(n,)+a2(n,,) • ^^-^^^ 

Npix is the number of pixels after the Galaxy cut and cr{h) is the r.m.s. error associated 
with the pixels in the direction h. The object hij is the reflection of hj in the plane with 
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normal fjj, 

hij = fij — 2{hi ■ hj)hi. (5.62) 

Lower values of S{hi) corresponds to a higher degree of symmetry. The temperature fluctu- 
ation is computed using the eigenmode expansion and a flat spectrum. The resultant 6T/T 
depends on six parameters: the three spatial orientations fixing the domain orientation and 
the three topology scales. 

The smaller a given direction the more extreme the asymmetry. In for instance, if 
the z direction is extremely small, then the size of fluctuations will be notably smaller in 
this direction leading to a map where there is essentially less and less structure at OOBE— 
DMR resolution in this one direction and the structure will be drawn out on the (x, y) plane. 
Similarly for there will essentially only be structure along the large direction and so the 
pattern will appear to be rings of fluctuations along the small directions. 

Using this statistic, and models were constrained to be greater than half the 
radius of the SLS in their small dimensions. Subsequently, the circle method was used to 
study specific asymmetric models where again pessimistic conclusions were reached as 



discussed in S 6.2.1 
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6 Observing hyperbolic topologies in the CMB 



The eigenmode decomposition is weh known on simply connected and is given in §4.2 
However, when hyperbohc space is fully compact the decomposition becomes intractable. 
This is a stronger statement than simply saying the eigenmodes and eigenvalues are difficult 
to find. It is actually formally impossible to write down the eigenmodes analytically. The 



boundary conditions are so intricate they resist decomposition |27]. 

The absence of an analytic solution to the eigenmode spectrum can be directly related to 
the incipient chaos on compact manifolds. On simply connected geodesies show the first 
critical ingredient for the onset of chaos; that is, extreme sensitivity to initial conditions. 
The geodesic deviation equation shows that two nearby trajectories diverge away from each 
other exponentially quickly: 

^ = -Kp^^^C^u^ (6-63) 

where = dx^/ds and C,^ is the separation of neighboring geodesies. Normal to the 
geodesic flow this becomes ||8^ 

'^^^^ = -MeNCN^\ (6.64) 

where k is the curvature. If k < 0, these geodesies diverge exponentially and a coordinate 
invariant Lyupanov exponent can be interpreted as A = Still, there is no chaotic 

motion on infinite since there is no mixing and no folding of trajectories. In other words, 
there is no loss of predictability as a result of the exponential deviation. By contrast, when 
the space is made topologically compact, the mixing and folding of trajectories is assured 
and the flows are well known to be fully chaotic. An entropic measure of the chaotic flow 
can be related to the volume of the manifold through the Kolmogorov entropy Sk oc V~^/'^ 



where V is the volume of the spacetime [B7|. Notice that if the space is infinite this entropy 
vanishes. 

Chaotic flows on compact 2-dimensional manifolds in particular have been studied at 



great length |27, 88, 9C, 91|. The cosmological implications of chaos have only been touched 
upon [46, 59 1 and remain a largely unexplored terrain. By and large, people have tried to 
obviate the chaotic flows entirely when analyzing the CMB. Methods include brute force 
numerical determination of the modes, the method of images construction of CMB maps, 
and geometric methods. We will discuss a catalog of such studies. 

6.1 Direct methods in hyperbolic space 
6.1.1 Simply connected hyperbolic space 

Since all of the numerical methods will rely on the expansion of the eigenmodes on the 



universal cover, we summarize those results here. Prom the standard expansion (4.25). We 



need to determine each of the factors in the decomposition. In negatively curved space 



perturbations are time dependent according to eqn ( 4.24 ) with solution 

5 (sinh^ Tj — 3r] sinh r/ + 4 cosh rj -\- A 
(cosh rj — 1 



F{V) = — — :t3 -■ (6-65) 
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In a negatively curved space the Helmholtz eqn becomes 



1 



sinh^ r 



dr (sinh 



+ 



sm 



-de {sin Ode) + 



1 



sm 



(6.66) 



As emphasized in Ref. [92|, a complete orthonormal basis for the simply connected space is 
formed by modes with real > 1 and so the eigenvalue range is k"^ = [l,oo]. These modes 
vary on a scale below the curvature radius and are therefore subcurvature modes. Standard 
causal theories such as inflation will not seed fluctuations on scales k < 1 although some 
other unforeseen mechanism may. It is customary to introduce another parameter 

g2 = - 1 (6.67) 
with the range = [0, oo]. The eigenmodes on can be expressed as 

tpqlm = Xgi{r)Yimin) (6.68) 
where the l^m's are the usual spherical harmonics and the Xg£ are the radial functions 

iq) 



X, 



qi 



r(£ + i 

T{iq) 



1 



sinh r ''•i 



P. 1^ (coshr) 



(6.69) 



where F denotes the Gamma function and P denotes the Legendre functions [p2| . 

The are drawn from a Gaussian probability distributions with proper normalization 
for hyperbolic space 



q{q^ + 1) 

We can summarize the Sachs- Wolfe effect as 

5T{h) 



V<s>{q)6{q - q')5u'5„ 



d''k^An)Lr 



with both the surface Sachs- Wolfe and the integrated Sachs- Wolfe accounted for in 



The Fourier multipoles are 

dk 



1 ^ 



+ 2 



drjF' {rf) 



Vsls 



Cf = 2^" 



fc(A;2 _ i; 



V^{k) 



1 Q no , 

-Xh(A7?) + - / dr Xk,{r)F'{Ar, - r) 



(6.70) 



(6.71) 



(6.72) 



(6.73) 



The COBE-DMRdata alone tends to favor 17 ~ 0.3 — 0.4 (see for instance [93|) although both 
very low O and an Q ~ 1 are still compatible. For a long time, the COBE results along 
with other astronomical observations put low cosmologies in the limelight. However, 
more recent high i observations [^] in conjunction with the supernovae data |64] have 
pushed public opinion back towards 17 = 1 cosmologies. The recent high i observations 
only examine small patches of the sky. A full sky map is required to advance topology 
observations and so we reserve more commentary until the launch and hopeful success of 
the future MAP and Planck missions. 
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Figure 31: Embedding of the cusp topology with one of the dimensions suppressed. 



6.1.2 Cusps 



There are multiconnected hyperbohc spaces which are not completely compact (see Ref. 
for some examples such as those of Refs. Q). For some of these the motion is not 
chaotic and it is possible to find the eigenmodes. One topology of particular interest is the 
toroidal horn p^, Bq]. The horn is interesting since many manifolds have cusped corners 



as stressed by |13]. Despite their frequency in a set of generic manifolds, it was argued in 
Ref. that according to the thick-thin decomposition, it would very improbable for us 
to live in the thin part of the manifold. However, this argument is flawed. Since the cusp 
narrows exponentially quickly, an observer can live in a fat part of the manifold and still see 
photons coming from a constricted part of the cusp. However, Olson and Starkman have 
developed an approximation scheme to study cusps on full compact manifolds which shows 
the generic patterns do in fact persist. 

The cusp is most easily understood in the upper half space representation of H'^ (see 
Appendix ^) 



-drf + dz^ + e 



2z 



{dx^ + dy' 



(6.74) 



with the coordinate transformation = cosh r— sinh r cos 0, e~^x = sin 9 cos (j) sinh r,e~^y = 
sin 9 sin (p sinh r . The (x, y) subspace looks like a conformally stretched flat space. The topo- 
logical identifications x — > x + L^ and y y + L+y make a 2-torus which is squeezed in area 
by the exponential factor e~^^ into a tightening cusp. Since the subspace is conformally flat, 
there are no tangled geodesies and it is possible to decompose the eigenvalues analytically. 



In the coordinate system (6.74), the temperature fluctuations can be decomposed as 



— (nj^)=y dq ^g^-'^y(^)^g^-^y^g^-^v 



(6.75) 



with the normalization 



N, 



2gsinh(7rg) (2 - (5„^o)(2 - Sr, 



1/2 



bh 



(6.76) 



and 



L. 



qn^ny 
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- + 2 
3 



Vo 



drjF' [t]) 




(6.77) 
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Figure 32: Flat spot in a map of the sky in the cusp topology with scales b = h = 1 and 
no = 0.3 



where Kiq is a modified Bessel function with imaginary index. The entire function is 
included here in the integration over r]. The argument of the Bessel function is = 
47r^(n^/L^ + Uy/Ly). Notice that q is still a continuous index since the z direction is 
infinite. There is no cutoff in the range of q but there is a cutoff in the range of (uxjUy). 
This cutoff results in a severe suppression of power as an observer lives nearer the cusp. 
The suppression appears as a flat spot in simulations of the COBE— dmr sky with concentric 
rings of larger and larger structures as the cusp widens |98]. 



A standard likelihood analysis of the C^s places the size of the torus at the location of 
the Earth to be > Arj. This is not a particularly narrow part of the cusp relative to the 
SLS but nonetheless an observer can still see exponentially deep down the throat putting 
large quiet regions in the CMB maps [p^j. 

The generic property of flat spots was tested in Ref. |9£]. They showed that from 



a typical location in a cusped manifold m003, flat spots of angular size ~ 5° would be 
visible for = 0.3. They conjecture that observable flat spots on this scale are typical of 
cusped manifolds. To handle the fully compact case Olson and Starkman found a means 
to approximate the modes in the cusp without tackling the full eigenmode solutions. They 
consider a horosphere, a sphere within the Poincare representation of H^, which is tangent 
to the cusp and passes through the point on the SLS. On the horosphere, the transformation 
group simplifies and they are able to isolate modes. 

The specific space they consider is a cusped manifold, namely m003, constructed from 
two tetrahedra with total volume V ~ 2.0299. The restriction to the horosphere results 



in a hexagonal tiling of and the eigenmodes are reminiscent of those of eqn ( 5.58| ). 
Consequently the minimum eigenmode from the hexagonal tiling provides an estimate for 
the extent of the flat spot. Situating the observer at a symmetric point in the manifold they 
find that fiat spots in the CMB from cusped regions subtend an average angle of roughly 
5° for = 0.3. Larger values of 0,q will naturally diminish the angular scale. They argue 
that in general cusped manifolds will show fiat spots as characteristic features in the CMB 
maps. 

6.1.3 Numerical Eigenvalue Solutions 

The eigenmodes can always be obtained by brute force numerically. Three groups have 
developed numerical methods to isolate a list of eigenmodes and eigenvalues for a small 



collection of manifolds |10C, 101, 102, 103, 104]. A comparison of the simulated CMB sky 
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Figure 33: Dirichlet domain for the Thurston space 



with the COBE— DMRdata shows the specific compact hyperbohc manifolds studied were 
consistent with the data. However ah of the statistical analyses compare only the C^'s and 
not the full correlation functions. Using the method of images to simulate the CMB, a full 
statistical analysis is performed by Bond, Pogosyan, and Souradeep as described in section 



§6.1.4 where more negative conclusions are obtained; namely that the spaces studied were 
inconsistent with COBE except for maybe one orientation of the manifold. It would be 
interesting if numerical eigenmode constructions would compare a full statistical analysis 
to that performed using the method of images to check for consistency. 

Eigenmodes of the Thurston Space 
Inoue was the first to find precise eigenmodes on a compact hyperbolic 3-space for the 
purposes of simulating a cosmological model. (Aurich and Marklof were the first to compute 
the eigenmodes on an orbifold jl05| ].) His method was based on the direct boundary element 
method initially developed by Aurich and Steiner for the study of 2-dimensional compact 
hyperbolic spaces. The numerical method is based on solving the Helmholtz equation 



(V2 + fc2)^,(x) = 



(6.78) 



on a compact manifold, i.e. with periodic boundary conditions on the universal covering 
space. Various methods for doing so include the finite element method and the finite dif- 
ference method. He uses the more precise but numerically more time consuming method 
called the direct boundary element method (DBEM) also used by Aurich and Steiner to 
study quantum chaos on a 2D compact space ||8^ . He isolated the first 36 eigenmodes of the 
manifold m003(— 2, 3) from the SnapPea census, otherwise known as the Thurston space in 
homage to the Fields Medalist W. Thurston. There are two particularly interesting results. 
Firstly he finds a cutoff in the low k eigenvalues and therefore a maximum wavelength at 
a value near the average diameter of the space. Secondly, he discovered that the expansion 
coefficients of the eigenmodes {iqim below) are well described as pseudo-random in behav- 
ior I IOC]. Despite the long- wavelength cutoff, the integrated Sachs- Wolfe effect is able to 
compensate for low Q.q, as expected. Consequently, assuming an Harrison-Zeldovich power 
spectrum, he finds that the C^s are consistent with COBE-DMRfor 0.1 < ilo ^ 0.6. 
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As in flat space, the eigenvalues are a discrete subset of the continuous eienvalues. Since, 
the eigenmodes of a compact manifold are a subset of the eigenmodes on the universal cover 
(§ |6.1.1| ), they can be expanded as 

V'^ = E ^gemXge{r)Yem{n). (6.79) 
The can be thought of us as containing the analogue of the relations found explicitly in 



flat space in section § 5.1.2. The challenge of numerically isolating the modes can be reduced 
to the still difficult task of finding the (,qem- Inoue explicitly finds that the S,qim behave as 
random Gaussian numbers for {£ < 19 and q > 9.94). He extrapolates that this continues 
to be true, even more so, for all higher modes and introduces this as an approximation 
to determine modes above k = W. A distinction should be made between the random 
Gaussian initial conditions contained in the ^q- and the random Gaussian behavior of the 
^g-. The former is an assumption which may or may not be true depending on the history 
of the universe while the latter is a property of the manifold. 

The DBEM was first used to find the pseudo-Gaussian random numbers £,q£m for A; < 10 



I IOC] and the analysis was then extended to include the first 36 eigenmodes with k < 



13 [101]. (As discussed below, it is claimed in Ref. |103] that 2 modes were originally 
overlooked in this method). A numerical cutoff in the low k spectrum was found at the 
value ki = 5.41. The corresponding maximum wavelength Amax = 27r/A;i can be compared 
to an approximation for the diameter of M which he computes as the average of the in- 
radius and the out-radius. The manifold has in-radius Tin = 0.535, out-radius rQ^f = 0.7485 
(and volume V = 0.98139). This yields k^i^ = 4:7r/{r^ + r_) = 4.9, an agreement to within 
10% of his ki. There are no supercurvature modes in the Thurston space. 

The pseudo-Gaussian behavior can be interpreted in a quantum mechanical context 
]10C]. Let the Laplace-Beltrami operator be the Hamiltonian of a quantum system, then 



the eigenmode is a wavefunction eigenstate. Since the underlying classical dynamics is 
chaotic, it has been conjectured that the quantum mechanical system will be governed by 
the predictions of random matrix theory (RMT) |106| , |9l[ ]. One such prediction is that the 
square expansion coefficients are Gaussian distributed if expanded with respect to a generic 
bases. In particular, the coefficients should be given by the Gaussian Orthogonal Ensemble 
(GOE) 

P{x) = -^e-^/\ (6.80) 
\/2'irx 

To check this prediction for the modes ( 6.78| ), Inoue first recasts the eigenmodes in terms 



of real independent coefficients aqim and real functions Rqim as 

= E ^qimRqem (6.81) 

£m 

where the a„^m can be expressed in terms of the and the Rqim can be expressed in 



terms of the Xq^Yi^ |10C| ]. Following ]88], he examines the statistical behavior of 

X = '"^'"^ 7 (6.82) 
^1 
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where aq is the average and the variance. The probabihty is singular at x = and so it 
is customary to compare the numerics to the cumulative distributions 

J(x) = r P{x)dx = erf(i/^). (6.83) 

To test the RMT prediction, the numerically determined cumulative distribution are com- 
pared for a goodness of fit. Very good agreement with the GOE prediction is found confirm- 
ing Gaussian behavior for the low lying states. The Gaussian behavior is expected for the 
highly excited states as a consequence of the classical chaos but Gaussian behavior for the 
low lying states is less obvious. Yet this is in fact what his observations confirm. He also 
tests the randomness of the a^^m and finds that they do behave as random variables. The 
randomness is not a property of the eigenmodes but is due rather to the almost random 
distribution of images in the universal cover. The number of copies of the fundamental 
domain inside a sphere with radius r/o is 

n, = ^(^^^M2r?o) - 2r,o ^^ ^^^ 

giving ni ~ 29 for r]Q = 1.6. 

The number of eigenmodes increases as and the number of boundary elements in- 
creases as k"^ . The task of computing high k modes becomes unmanageable with this 
method. For higher A:, Ref [|100(| suggests taking the clue from the behavior of low k modes 
and approximating the highly-excited states as random Gaussian numbers with a variance 
proportional to q^"^. Using Weyl's asymptotic formula, 

N[q] = -\ q = V/c^ - 1 q » 1 (6.85) 
67r'^ 

with an integer. Since for a CH space the volume is fixed, Weyl's formula effectively 
relates the number of states at a given k with a topological feature of the space. It also allows 
an estimate of the kj. The spacing between discrete eigenvalues decreases as the inverse 
of k for large k and so approaches a continuous spectrum in the large limit, as expected. 
Weyl's formula is well obeyed and 30 of the 36 modes show random Gaussian behavior. 
Interestingly, six degenerate states are found which correspond to a nearly symmetric mode 
reflecting the global symmetry of the fundamental domain. While a linear combination of 
the degenerate modes shows Gaussian behavior again. This is typical of classically chaotic 
systems where the classical chaos leads to a Gaussian behavior in the quantized eigenmodes 
although occasionally the global symmetry of the space can surface in the eigenmodes as 
non-Gaussian behavior. (The connection with quantum chaos led Inoue to use the Selberg 
Trace Formula to compute eigenvalues for a large number of CH manifolds |^] . This method 
proved to be quicker and easier to implement numerically.) 

The C^'s are calculated by sewing together the numerically obtained eigenmodes for the 
eigenvalues in the range 5.4 < A; < 13 with the above approximation for the expansion 
coefficients for 13 < A; < 20. With the assumption that =constant, the resultant is 
compared with the COBE-DMRdata to conclude that the spectrum is consistent with COBE— 
DMRfor 0.2 <VtQ< 0.6. The cutoff in the spectrum due to the minimum mode k = 5.4 is 
buried under the contributions to low d. power from the integrated Sachs- Wolfe effect. The 
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ISW becomes important during the curvature dominated epoch at about 1+z ~ (1— $7o)/f^o- 
They find the contribution from A; < 13 modes to Ci for 2 < £ < 20 is 7 % for 0, = 0.2 and 
10% for Qo = 0.4. The rest is due to the ISW. For this reason, the spectra appear to have a 
gradual peak near the long-wavelength cutoff or are nearly flat down to low k and are less 
severely constrained for Q > 0.1 than the flat models of section §|5| [100, 101]. In fact they 



find very good agreement with COBE-DMRfor JIq = 0.6, better than other FRW models. In 
particular, an alignment of the peak in the COBE— DMRdata at ^ ~ 4 with the peak due to 
the cutoff accounts for the better fit. 

The conclusion cannot be stated as proving that the Thurston space is consistent with 
the observations, only that it does not appear inconsistent. A comparison of the full corre- 
lation function to the data may not survive consistency and in fact does not according the 



method of images analysis of |17, 18|. 



Eigenmodes of an Orbifold 

Aurich was able to compute a huge number of eigenmodes, the first 749, in an orbifold 
with negative curvature. Orbifolds can have compact volume but possess points which are 
not locally R^. They can also have rotation group elements among their isometrics, unlike 
manifolds which have no groups with fixed points in their isometrics. In the absence of 
a predictive theory for the topology of the cosmos, there is no convincing reason to omit 
them from the catalog. Interestingly, they can have even smaller volumes than the minimum 
bound on compact hyperbolic spaces. 

The fundamental cell is a pentahedron which is symmetric along a plane dividing the 
fundamental domain into two identical tetrahedra. This allows a desymmetrizing of the 
pentahedron useful for deducing the eigenmodes following early work |105[| . There are 9 
compact hyperbolic tetrahedra and the one built into the pentahedral orbifold is known 
as Tg. The volume of Tg is Vtetrahodron — 0.3586524 which is smaller than the Weeks 
space. The smallest compact hyperbolic tetrahedron T3 has volume V = 0.03588506, ten 
times smaller than Tg and smaller even than the existing bound on CH manifolds {Vjn = 
5/(2\/3)arcsinh^(\/3/5) ~ 0.167). Using the boundary element method 749 eigenmodes on 
Tg are obtained. 

Both radiation and matter are included in his analysis which evolves the metric per- 
turbations according to the usual adiabatic, linear perturbation theory and solves F{r]) 
numerically. In his analysis, the eigenmodes are expanded as 

$(77,x) = 5^F,(r?)V^,(x) (6.86) 

where the ipg are the eigenmodes and all of the time dependence is in Fq{r]). Notice that in 



comparison to the standard expansion 4.2, he does not include the usual randomly seeded 



fluctuation amplitude ^t- Instead all of the initial conditions are absorbed into F{r]) with 

Fq{m) = Ki'li) = 0. (6.87) 

There is however some element of randomness in the sign of the fluctuation. The eigenfunc- 
tions have a freedom in the phase. Since the phase is always chosen to be real, this freedom 
results in a lb ambiguity in the eigenfunction which is chosen randomly. The constant a is 
normalized against the COBE-dmr data. This is consistent with the Harrison-Zeldovich flat 
spectrum in terms of normalization but lacks any of the randomness that real fluctuations 
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Figure 34: This figure was supplied courtesy of R. Auricli. Tlie figure shows measured 
in micro-Kelvin versus ^ as computed by Aurich and Steiner for an JIq = 0.9 cosmology 
with r^A = 0.6. The angular power spectrum Ci for the orbifold is shown in comparison 
to an infinite topology as obtained by CMBFAST. The suppression for ^ > 30 is only due 
to the truncation of the sum over the eigenmodes since only the first 749 eigenmodes have 
been taken into account. 



would have. Because the randomness of the temperature fluctuations about the mean is 
neglected, it would not be quite right to interpret the simulated maps of Ref [104] in terms 
of an actual universe. The advantage however is that the maps show a random nature which 
can only be due to the random character of the modes themselves. This is important for 
the reasons discussed at length in the previous subsection |6.l| . 

Since the produces a global averaging, an ensemble average is generated by the expan- 
sion (|6.86| )- (|6.87 ) even if the standard random initial fluctuations were ignored. Although 
the simulations depend on the location of the observer and the orientation of the manifold, 
this dependency is not examined and the observer is fixed, offset from the origin. The C^'s 
are compared to the £ < 30 COBE— DMRdata, the Saskatoon data around I ~ 100 | 107 | and 
the QMAP data [108] above I ~ 80. All modes are computed up to kmax = 55. Aurich 
varies f^o between 0.2 and 0.6 and finds reasonable agreement for JIq — 0.3 — 0.4. For 
smaller VLq he finds that the numerical Ci increases too quickly with i. We suggest that this 
is evidence of a long wavelength cutoff in the spectrum. Although for a lower Oq, the ISW 
is important in terms of building up the low d. power, it is also true that the suppression 
is more severe. To put it another way, if the constant a were normalized to COBE-DMRat 
high ^, as it should be in order to minimize the effects of cosmic variance, then the steep 
slope would be seen as a suppression of large angle power due to the finite extent of the 
space. 

The same orbifold was studied again in Ref. [109] where both a cosmological constant and 
a smooth dark energy component were added to the radiation and matter energy density. 
The total energy is taken to be nearly but not quite flat to address the recent results from 
the small-angle CMB experiments ]|38|, pq]. They still find a suppression in power for i < 10, 



as illustrated in figure 34 



Although the results of Ref. 1104] to be at variance with the results of [0, ^ we 
emphasize again the limitations of comparing only the C^'s. This is expanded upon in 
section §6J^. 

Fast Method for Isolating Eigenmodes 
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Figure 35: Dirichlet domain for the Weeks space 



A fast numerical method for obtaining the eigenvalues on compact hyperbolic manifolds 
was developed by Cornish and Spergel |103| ]. In comparison to the boundary element method 
developed by Aurich and Steiner ||8^, it is technically inferior but far quicker and easier to 
implement numerically. The DBEM has only been applied to a few cosmological spaces, 

and the Weeks manifold to 
allows them to obtain the 



the tetrahedral orbifold of [104|, the Thurston manifold jlO' 



name a few |8£], to name a few, while the method of Ref. |10£ 
lowest eigenvalues and eigenmodes for 12 manifolds, requiring only the generators as input. 
The method is to simply solve 



using a singular value decomposition. The eigenmodes are again expanded in terms of the 
eigenmodes on the universal cover as in eqn ( |6.79| ). Random points are selected within 
the fundamental domain and all of the images out to some distance are located in the 
covering space. For each point pj there are nj such images located with the generators ga- 
The nj{nj + l)/2 boundary conditions are 



^k{gaPj) = "^kigpPj) ai^ (3. 



(6.89) 



Using the expansion of the eigenmodes in terms of the eigenmodes of the universal cover 
|6.1.l| , they write (|6.89| ) as a collection of difference equations which they then solve using 
a standard singular value decomposition method for handling over constrained systems of 
equations. 

The example they studied in detail is the Weeks space m003(3,-l). They also looked at 
the Thurston space and found agreement with the DBEM method applied by Inoue except 
they found 2 new modes missed previously but later confirmed (6^]. Again, for higher /c, 
the number of modes obeyed Weyl's formula well. Also, the modes were well described by 
random matrix theory which predicts the expansion coefficients obey a Gaussian Orthogonal 
Ensemble. 

The Weeks space has a large symmetry group which ensures a large degeneracy of 
eigenmodes. The symmetry group is the Dihedral group of order 6. They found the first 74 
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Figure 36: This figure was supplied courtesy of N. Cornish. The lowest eigenmode in the 
Weeks space is shown. The space is drawn in the Poincare representation. The three 
panels represent different slices through the Poincare ball. The upper figure shows the 
slice through the fundamental domain while the lower figure shows the slice through the 
eigenmode. Specifically, the leftmost panel is the x = slice, the middle is the y = slice, 
and the rightmost panel is the z = slice. 



eigenmodes, many of which are degenerate. The higher the mode, the more degeneracies in 
agreement with Weyl's asymptotic formula. A view of the lowest eigenmode in the Weeks 
space can be found in figure 

Implementing their eigenmodes in a simulation of a cosmological model, they compare 



the numerically generated C^'s to the COBE-DMRdata |103]. Their numerically modeled cos- 
mologies are based on the first 100+ modes of the following small spaces: The Weeks space 
m003(3,-l) with V = 0.9427, the Thurston space m003(-2,3) with V = 0.9814, s718(l,l) 
with V = 2.2726 and v3509(4,3) with V = 6.2392. One realization of the Weeks space 



is shown in figure 37. Generically, a kind of statistical isotropy prevails even though the 
spaces are globally anisotropic. Since the expansion coefficients are pseudo-random, Gaus- 
sian distributed numbers which are statistically independent of I and m, they generate 
nearly isotropic eigenmodes and lead to a kind of isotropy across the microwave sky. The 
authors suggest inflation may not be need to explain why the universe is nearly isotropic. 
However inflation is still needed to explain why local values are so marginally near flat. 

A standard likelihood analysis based on 100 realizations of each topology was performed. 
A mode cutoff in the spectrum is again evident and gets worse as JIq is lowered. For all of 
the manifolds studied, they find that 

= (1.3 ^ l.6)D (6.90) 

where = 2'ir/qi (although it is customary to define the wavelength corresponding to a 
given mode as A'^' = 27r/k) which obeys the bound 

4-D ,0 . /27r^2 



L>2(sinhL> + L>)2 \DJ 



<k'i<l+[^) (6.91) 
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Figure 37: This figure was supplied courtesy of N. Cornish. A numerical realization of the 
CMB in the Weeks topology with f^o = 0.3. One pair of matched circles is indicated by the 
white lines in this projection. 



with D the square root of the smallest integer > D^. They do caution that they cannot 
find modes with q = [0,1/4] using their method, and so perhaps supercurvature modes 
lurk. However the ISW becomes increasingly more important as r^o is lowered. There is 
an optimal value where the two effects compete to create a spectrum with a slight tilt at 
low i and give a better match to the COBE— DMRdata. For the Weeks' space this occurs at 
= 0.3. 

Additionally, once the spectrum is normalized to COBE— DMRthey find that the size of 
fluctuations on 8h~^ Mpc is naturally increased above the value of ag = 0.6 for a simply 
connected universe with = 0.3 to a value of erg = 0.75 for the Weeks manifold with 
= 0.3. The present day cluster abundance seems to imply = 0.9 ± 0.1 for JIq = 0.3 
|11C|] and so higher values are desirable. 

They find the C^'s have a better fit to the COBE-DMRdata with a relative likelihood of 
~ 20 and a better fit to the large-scale structure date (o"8 increases by ~ 25%). However, 
as with all of the other statistical analyses they have applied a weak test by analyzing 
the angular averaged Cis and not the full correlation function. Their conclusions can be 
interpreted as finding the manifolds are not ruled out by the Ci alone. 



6.1.4 The Method of Images 

The fluctuations can be simulated using the method of images. The correlation functions 
can be calculated without the explicit eigenmodes and eigenvalues. The procedure is to 
sum the correlation function on the universal cover over all images out to some large radius. 
The more distant images are handled in a continuous approximation. The eigenspectrum 
can be calculated but the correlation functions are obtained to better accuracy for a given 
order in the sum. Bond, Pogosyan and Souradeep implemented a detailed method of images 
|l7|, fsl. 111, 112, 113|. They emphasize that the full correlation function C{n,h') must be 



compared to the data for a meaningful statistical analysis. The philosophy is to perform 
a statistical search for patterns. They find two principle effects (1) anisotropic patterns 
and (2) a long wavelength cutoff in the power spectrum. The patterns appear as spikes of 
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positive correlation when a point on the SLS and one of its images is correlated. The larger 
the SLS relative to the out-radius, the more statistically significant will the patterns be. If 
the SLS is smaller than the in-radius then naturally the correlation function is very close 
to that for the simply connected space. 

The angular correlation function C{fi,h') is computed from the spatial two-point corre- 
lation function = {^{x,tls)^{x' ,tls)) which can be expressed as 



i j=l 



(6.92) 



with the sum over a discrete ordered set with multiplicities mj. Notice that the obey 

(6.93) 



as always, however with the slightly different two-index notation. The j index accounts 
for any degeneracies. The spatial correlation function on the compact manifold M'^, 
can be expressed in terms of the spatial correlation function on the universal cover. A 
derivation of the relation between ^'^ and exploits orthonormality and completeness with 
the following equations: 



/ dx'C$^(x,x')*^,(x') =P$(fc,)*^,(x') 
Jm 

[ dx'e<i."(x,x')^'"(A;,x') =p4.(A:i)^"(A:,x'). 



(6.94) 



(6.95) 



Since eigenfunction on the compact space must also be eigenfunctions on the universal cover 
(although the converse is not true) it follows that 



/ (ix'e$"(x,x')^-(A:,x') 
JM" 



j dx'C4.^(x,x')^^,(x'). 



Combining (3.96) with ( 6.95| ) gives 

/ dxV(x,x')^5(A:,x')= / dx'C$"(x,x')^?,(^') 
JM JM^ 



(6.96) 



(6.97) 



Since M. tessellates M."^ we can re-express the left-hand-side of (|6.97|) as 



I dx'e$"(5,5x')n(^')= / ^x 
^gpjx Jm 



E?<f"(^'9x') 



*f,(x') 



(6.98) 



where in the last step a regularization is needed and denoted by the tilde above the sum- 
mation. Identifying integrands gives 



^$^(x,x') = 5]e<l."(i,5x') 

96r 



(6.99) 
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and the ^^'^ can be calculated as the sum over images in the universal cover with only a 
knowledge of the group elements. The spatial correlation function on the universal cover is 
known to be 

- = r T^'-^V^{<D (6.100) 

Jo (g^ + l)gsmhr 

and is taken to be a Harrison-Zeldovich spectrum. This is the core of the method of 
images. 

The regularizer requires some additional effort. The correlation function on the universal 
cover does not have compact support: 

/ (ix'^$''(x,x') = cx). (6.101) 

The need for regularization is not a result of the CH space having a large number of periodic 
orbits as incorrectly claimed in nor is it dues to the chaotic nature of the trajectories but 
rather is a result of the correlation function not having compact support. Even flat spaces 
require regularization. The regularized spatial correlation function can be written 

iKx, x') ^ e<l."(x, X ) - / (ix"e$"(x', x") (6.102) 

vM JgM 

for g such that x' lies in gM. The regularization ensures 

dx|^(x,x) = 0. (6.103) 
Finally then, 

e*'=(x,x') = ^ei(x,5') = ^e$"(x,x') - ^ / dx"C*"(x',x"). (6.104) 

The regularization prescription is not unique. The actual limiting procedure utilized in 
Ref. fl^ sums images up to a radius and then regularizes by subtracting the integral of 
^$"(r) over a spherical ball of radius r^,: 



^$'^(x,x') = lim 



^ C^^ir,)-— / drsinh^ re$"(r 

r,<r. Jo 



(6.105) 



with rj = d{x,gjx') < ^j+i. The value of r^, is numerically pushed out to 4 or 5 times the 
out-radius in order to get a convergent result. This procedure is simpler and does well for 
large r,,. It also does not require a detailed dependence on the complicated shape of the 
fundamental domain. 

An example of the maps they generate using the method of images is shown in figure 



38 



They are also able to estimate the density of states and therefore obtain a rough estimate 
of the power spectrum using the method of images. This is in addition to their primary 
result of having computed the correlation function. When they do estimate the density of 
states they always find an infrared cutoff as expected for long wavelength modes. 
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NOP PGP 




Figure 38: This figure was supplied courtesy of T. Souradeep. Tlie full CMB sky is rep- 
resented by the two hemispherical caps - one in the direction of the South Galactic Pole 
(SGP) and the other in the direction of the North Galactic Pole (NGP). The label SHC 
refers to a small compact hyperbolic model, namely m004(-5,l). As well as showing the 
fluctuations, this figure shows correlated circle pairs explained in more detain in a following 
section. 



This is consistent with the expectations. For a 3-dimensional CH space there is a bound 



on Cheeger's isoperimetric constant of |114, 115 ] 

1 



krain > hc/2 > 



1-1/2 

2 / dt cosh^it) 
Jo 



0.92/dM (6.106) 



and so there are no supercurvature modes for dju < 0.92. The suppression of power is 
covered in part by the ISW and so is less prominent than in flat models. 

Given their numerical calculation, they do a Bayseian probability analysis comparing a 
few compact hyperbolic cosmologies to the COBE— DMRdata. Their analysis represents by 
far the most complete of the statistical tests that has been performed and in principle is 
the most complete test that can be performed. They find the compact hyperbolic models 
they study are inconsistent with the COBE— DMRdata for most orientations although for a 
small set of special orientations they do find a better fit than the standard infinite models. 
Their results are qualitatively independent of the matter content. The analysis magnifies 
the inadequacy of using the Ci alone. In particular they point to the large cosmic variance 
that results from the break down associated with isotropy. 

With the assumption of Gaussianity both in the noise and in the raw 5T/T, the proba- 
bilistic comparison used is 

where Ct = C{fi,n') and 5 is the data. To numerically obtain the full correlation function 
requires an evaluation of the spatial correlation function at Np(Np + l)Nl/2 pairs of points 
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where Np is the number of pixels and Nl is is the number of points along the line of 
sight used to integrate the ISW. They find an Nl ~ 10 is satisfactory. They estimate the 
likehhood function 

£{Ct) = r{6\CT) = J d6r{6\6)r{6\CT) (6.108) 

with 5 the map which maximizes the conditional probability and the integration is performed 
over all realizations of the simulated sky 6. The resultant likelihood function they use is 

(2^)^./2||C7^ + C7^||i/2" ^^-'^^^ 

with Cat the noise covariance matrix. What is actually obtained is model dependent relative 
likelihoods. The model-dependent parameters are M, the orientation, the location of the 
observer, and ^}\, and the initial assumptions regarding the spectral shape and 

character. The manifolds studied were m004(-5,l), which is relatively small, and v3543(2,3), 
which is comparatively large. (JIa = always. Adding relaxes the constraints.) They 
varied = over three values arranged so that the SLS was comparable to the out- 
radius and varied over 24 different orientations. They do leave the observer at the maximum 
of the injectivity radius. This location gives the most symmetric perceived shape to the 
Dirichlet domain which one might expect to lead to the most conservative bounds. It seems 
fair to assume that moving the observer to a thinner region in the manifold for instance 
will only amplify asymmetries and the constraint on the correlation function should only 
be more severe. However, others have argued that locating the Earth away from the local 
maximum of the injectivity radius can increase the likelihood fit of a CH model to the data. 
For instance, to check the effect of the inhomogeneity, Inoue moved the observing point 
for the same model m004(-5,l)=Thurston manifold and performed the Baysean analysis as 
Bond, Pogosyan, and Souradeep did. He found that there are a few choices of the position 
and orientation for which the likelihood is much larger than that of the infinite counterpart 
|l62|| (as did |]l8|). The best-fit positions are scattered in the manifold and far from the local 
maximum of the injectivity radius. 

In any case, the results of [0, IS] are presented as the relative likelihood of a given model 
in comparison with a simply connected hyperbolic CDM model with the same = ^o- 
They uniformly find that the compact hyperbolic models have very small relative likelihoods 
as compared with the simply connected models. The rare exception occurred near Arj ~ 
r+ with a particular orientation. The statistical significance is unclear since a fortuitous 
alignment of the measured fluctuations with simulated topological images could enhance 
the likelihood when taken over all realizations which reinforce this correlation. The authors 
defer conclusions to future tests such as the circle method. Although again, Inoue argues the 
statistical significance is clear. If the C^s fit the data poorly, there will be no orientation of 
the manifold that leads to an alignment of measured fluctuations with topological images. 
Put another way, it may be unfair to suggest that a good fit to the data is a fortuitous 
alignment instead of acknowledging this good fit to be a good fit. 

Bond, Pogosyan, and Souradeep draw the general conclusion that dM/"^ > 0.7 Arj. The 
compact models they tested are excluded at the 3a" level with the exception of those with 
very special orientations. By contrast, if they were to only analyze the statistical likelihood 
of the Ci they would mistakenly conclude that the compact models were preferred at the la 
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level over the simply connected cosmology. They emphasize that the error bars on the C^'s 
are huge because of the exaggerated cosmic variance in the topologically connected models. 
It may also be worth noting here that there is some argument over the interpretation of 
these error bars (see Ref. p^). 

Due to the global breaking of homogeneity, they find the variance is spatially depen- 
dent. This is another reason to be weary of conclusions based on Ce alone. They found 
characteristic loud spots, that is regions in the sky with larger variance than others. Loud 
regions correspond to intersections of the SLS through smaller regions of the fundamental 
domain. The regions are geodesicaly small, that is to say, there are shorter geodesies relative 
to other regions in the volume. Unless the region is very very small, the ISW can obscure 
this particular feature. These loud and quiet regions are familiar from the cusp topology 
studied in Ref. [98|. There, even the ISW cannot compensate for the deepest regions of the 
cusp. 

Another attitude to take would be to argue that if we do live in a compact hyperbolic 
manifold, some orientation is necessarily going to be much more likely, namely the right 
one. The special status of the 'correct' orientation for our manifold, assuming it is compact 
hyperbolic, was taken seriously by one of the pioneers in cosmic topology, Helio Fagundes 
|116| , |117| . He reasoned that the particularly significant hot and cold spots in the COBE 
maps found by Cayon and Smoot |118[| may shed light on the orientation of the manifold. 
These loud spots may be patches of high and low regions in the physical density, as opposed 
to being just statistical fluctuations. Since the density fluctuations eventually evolve into 
large-scale structure, then the hot and cold spots should correspond to physical superclusters 
and voids respectively. Fagundes used this idea to try to match sources in catalogs of 
galaxy superclusters and voids and thereby flx the orientation and location of the earth in 
the manifold. His algorithm begins with a point P' centered on one of the spots isolated 
in Ref. [118| as illustrated in the schematic flgure 39. The point P' £ SLS. He then 
maps ghost images of P' = jP for 7 G F so the 7 are composite words built from the 
face-pairing matrices. He tries to alien the image points with voids and superclusters for 
different orientations of the space and different basepoints, i.e. different locations of the 
earth. He performed these scans for the ten smallest compact hyperbolic manifolds known 
III]- 



6.2 Geometric methods in hyperbolic space 

Given the likelihood analysis of the previous section, should we generically conclude that all 
compact hyperbolic spaces must be large relative to the observable universe? The answer 
is really "no" since there is so much model dependence in any of the direct methods of the 
previous section. We cannot be confldent about our assumptions for the initial fluctuation 
spectrum, the choice of manifold, the orientation, location of observer, local parameter 
values or even the statistics themselves. Remembering that there are an infinite number of 
manifolds to consider and the ambiguities in the model parameters, statistical conclusions 
are limited to the specific. These attributes beg for a template independent method of 
searching for topology. Much like gravitational lensing, we might hope to simply look at 
the sky and see evidence of topological lensing without prior assumptions about the shape 
of the lens. There is hope for such model independent observations as exemplified in the 
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Figure 39: This figure was supplied courtesy of H. Fagundes. It sliows the fundamental 
polyhedron as well as a copy of the fundamental polyhedron that intersects the spherical 
surface of last scattering. The point of intersection identifies the original location of a CMB 
spot in the universe. 



circles in the sky described below. Generically, we distinguish these geometric methods 
which search for patterns from statistical methods which rely on a specific model. 



6.2.1 Circles in the Sky 

Possibly the nicest geometrical observation made thus far has been the prediction of circles 
in the sky of Cornish, Spergel, and Starkman |2|. It has quickly become a popular 

topic in conversations on the topology of the universe. The circles are most easily seen in 
the tiling representation of a compact space. Each copy of the Earth will come complete 
with its own surface of last scatter. If the two images of the Earth are near enough, these 
identical copies of the SLS will intersect. The intersection of two spheres occurs along a 
circle. Since the observers at the center of the intersecting spheres are actually just copies 
of one observer, the circle of intersection must always come in pairs as illustrated in fig. 



40. To emphasize, we will not look up in the sky and see intrinsic circles in the microwave; 
that is, the circle pairs do not have identical temperatures along a circle but rather the 
temperature varies identically when taken along correlated circles M . An illustration of the 



location of circle pairs in a finite flat torus is shown in figure 41 from Ref. [120|. 

The most important aspect of the circles approach is that it applies to all multiconnected 
topologies and does not require a template nor any a priori assumptions about a model. 
All compact spaces will have circles if any part of the geometry is smaller than the SLS. 
The radius, number, and distribution of the circles will vary from space to space and hence 
the topology can be reconstructed from the circle pairs. A distribution of clone images in 
the Thurston space can be seen in figure Since the shape of the Dirichlet domain is not 
a topological invariant but rather depends on the location and orientation of the observer. 
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circles in the sliy 



the intersection of two spheres 
occur!! along a circle 




in the tiling picture 
one copy of the earth sees a circle to the lett 
and another copy sees a circle to the right 

from the earth we most see two identical circles 
one to tlic left, one to the right 



Figure 40: 




Figure 41: This figure was supplied courtesy of N. Cornish and shows the northern and 
southern hemispheres of the microwave sky in flat hypertorus. There are 13 matched circle 
pairs indicated by black lines. 
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Figure 42: This figure was supplied courtesy of N. Cornisli and illustrates just how many 
circles are to be expected. A fundamental domain for the Thurston space is shown. The 
distribution of points mark our clones out to a distance of three times the curvature radius. 
The large points are within one curvature radius, the medium sized points are within two 
curvature radii and the small points are within three curvature radii. Each of these clones 
will have a clone surface of last scatter. The clones which generate circle pairs will depending 
on the size of surface of last scatter compared to the clone distance and therefore depends 
on the value of Oq- 



even the Earth's location in the universe can be deduced. 



A statistical scan of the sky must be performed to draw the correlated circles out of 
the maps. The circle pairs will be completely hidden. To pull them out, consider two rings 
each with angular radius a and with relative phase cf)^ centered on arbitrary points x and 
y. To test whether these arbitrary rings are in fact correlated circles of intersection, the 
comparison statistic 

has been proposed where <>= J^^ d(p with S range [—1,1] [||. Perfectly matched circles 
have 5 = 1 while an ensemble of uncorrelated circles will have a mean value of 5 = 0. 



(Roukema includes the small scale Doppler effect in a slight alteration of the statistic [p5[.) 
Orientable topologies will have clockwise-anticlockwise correlations while non-orientable 
topologies will have a mixture of clockwise-clockwise and clockwise-anticlockwise correla- 
tions. In flat space, matched circle pairs have angular radius 



a = arccos 



with X the distance between the Earth and its image. In hyperbolic geometry 

( V (6.112) 

V sinh X tanh Arj J 



a = arccos 
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There are no circle pairs if the image is too far for the spheres to intersect and the expressions 
are invahd for X > Arj. The number of images grows exponentially with X in hyperbolic 
space so that most circles have small radii, although the statistic works best for large circles. 

Noise will degrade the circle statistic so that S ^ 1. The experimental noise in each 
pixel can be approximated as random Gaussian noise 

P(n) = — ^e-"'/2-^ (6.113) 

with variance Cn- The true temperature fluctuations are also taken to have a Gaussian 
distribution with variance 0"^. The probability for the matched circle is then 

with r(A^) the gamma function and ^ = (Ts/ct„ is the signal-to-noise of the detector while 
is the number of pixels [?]. For large N, the distribution has a maximum 

S^.. = Y^ + 0{N-'). (6.115) 

The higher the resolution and the signal-to-noise, the better this statistic will fare. If the 
experimental angular resolution is 69, then there are N ~ 2TTsma/56 data points around 
each circle of angular radius a. For COBE-dmr , the signal-to-noise ratio is = 2, 56 = 10° 
and N ~ 36 sin a pixels while for MAP, ^ ~ 15, 69 = 0.2° and N ~ 1800 sin a in its highest 
frequency channel. Using MAP parameters, only circles with a > 4° are detectable. 

To make a relative comparison with a probability distribution for unmatched circles 
they advocate 

P''iS)dS = ^^^^{1 - SY^'-'dS (6.116) 

which is centered at 5* = with a FWHM~ (81n2/n)^/^. The unmatched probability 
distribution was derived by assuming the temperature at each point is an independent 
Gaussian random variable. There are some weaknesses in this assumption since there are 
known correlations in the CMB. 

Cornish and Spergel applied the circles test to maps generated with their numerically 



isolated eigenmodes of section 6.1. For a realization of the Weeks space they found the un- 
correlated ISW clouded the statistic resulting in a poor match for circle pairs. To combat 
this pollution they remove all power in modes below i = 21 and find a substantial im- 
provement in the match with values of S" ~ 0.9. Their expectations for the future satellite 
missions are high. 

The circles method was first applied to simulated maps by Bond, Pogosyan and Souradeep 



using the method of images of section § |6.1.4| . From their correlation function they are able 
to make full sky maps. Doing so they search for circles using the cross-correlation coefficient 
between fluctuations along two circles Ci and C2 = gCi, 

Pu . ('^^^^9^^^ (6.117) 
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with X G Ci and ^rx G C2. A statistical average is implied by the angular brackets. For 
one given realization an integration over x along the circles is performed instead. They find 
P12 is in the range 0.6 - 0.95 for Qo = 0.9 and m004(-5,l) and 0.2 - 0.6 for S^o = 0.6 for 
v3543(2,3). Notice that is selected for each manifold so that V is comparable to the 
volume of the SLS. The matches they find are good even at COBE-dmr resolution. The 
correlations along circles gets worse as the ISW contribution is enhanced at low Qq. 

In addition to the circular intersections of the boundary of the SLS, there are inter- 
sections of the volume. Let S represent the collection of points contained within the SLS. 
Then the intersection of gS U S defines a lens-shaped region (see their figure) in the volume 
and must be identical to the intersection of S U g^^S. As long as r„ < At], then there 
will be lens-shaped intersections of the volumes and circular intersections of the copies of 
the SLS. Since the anisotropy occurs along the entire line of sight, Bond, Pogosyan and 
Souradeep consider the correlations in the full lens-shaped volume. Based on the COBE— 
DMR constraints obtained via the method of images as described in section §3.1.4 they took 
the topology scale comparable to the SLS. They argue that correlations must then be very 
near the faces of the Dirichlet domain and so the correlations provide a quick sketch of the 
shape of the universe |18]. 

The volume intersections are relevant since, contrary to intuition, they found the ISW 
is not entirely uncorrelated. While there are an infinite number of lines of sight which share 
at least one common point, there are also pairs of lines which have segments in common and 
so an enhanced correlation. Every pair of lines of sight which are directed toward the center 
of matched circles, necessarily contain segments of identical points (see fig.) and leads to 
correlated patterns. These substantial anti-correlated features tend to lie at the centers of 
matched circles. The anticorrelation comes from the interference term between the surface 
Sachs- Wolfe effect and the ISW. As described in section §3.1.4, their likelihood comparison 
of a handful of spaces to the existing data led to pessimistic conclusions. 

Another application of the circles method using COBE-dmr sought to identify an asym- 
metric fiat 3-space [p5| . COBE— dmrIs not ideal for detecting circles in hyperbolic space since 
the ISW accounts for most of the power in the COBE~dmr range of detection. However, 
for flat spaces with no cosmological constant the sky is determined by the surface Sachs- 
Wolfe effect alone without the obscuring effects of the ISW. Despite the low resolution, 
COBE-DMR might still see circles if we live in a compact fiat space. 

Roukema considered a specific asymmetric torus which was put forth as a candidate 
model to match cluster observations |121, 122]. Although he used the recent circles pre- 
diction to test his hypothesis, the symmetries of the space and the geometric approach 
are reminiscent of the earlier work of |^2| described in section §^]^. Hot X-ray bright gas 
in large galaxy clusters were used to search for topological images [121]. Two clusters at 
redshifts z ~ 0.4 very nearly form a right angle with the Coma cluster with very nearly 
equal arms. On the basis of this geometric relation, they take a toroidal geometry for the 
universe to explain these clusters as topological images of the Coma cluster. The distance 
from the Coma cluster to CL 09104-1-4109 and the distance from Coma to the cluster RX 
J1347. 5-1145 are both « 960/1-1 ]yjp(, fQj, 

zero cosmological constant. The third dimen- 
sion is taken to be larger than the diameter of the SLS and hence topological effects from 
this direction are essentially unobservable. The size of the small dimension was taken to 
be roughly Ary/13.2. Using the circles statistic, the candidate was ruled out at the 94 % 
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confidence level, provided that the cosmological constant is zero. Specifically, the statistic 
Roukema used was 



d^{- ' ..1/2 ) (6-118) 



A(f 



+ 



T 



where [A((5T/T)]j are the observational error estimates on the temperature fluctuations 
{5T/T)i as estimated by the COBE team. This statistic directly tests the consistency of 
temperature values within observational error bars |85]. If the ISW effect is treated as noise 



as in Ref . |85[| , the application of the circles principle using the statistic 6.118 can still enable 
rejection of a specific partly compact model or show that a model one tenth of the horizon 
diameter is consistent with the COBE data. 

In general, there are obstacles to implementing the circle method in practice. For one, 
the much emphasized ISW effect is not correlated in this way and so can obscure the circle 
pairs. Other problematic effects include the velocity and thickness of the SLS. In Ref. 
p5| , it was noted that many circles are partially lost along with 20° galactic cut and that 
both detector noise and foreground contamination posed difficulties for the circle detections. 
Realistically it may not be possible to observe circles even if they are there and this will be 
the challenge faced in realistic analyses of the future satellite data. 

6.2.2 Pattern Formation 

Developmental biology and condensed matter physics have long exploited pattern formation 
induced by periodic boundary conditions. Since compact topologies can be understood as 
a set of intricate boundary conditions, pattern formation has a cosmological analogue. The 
emergence of patterns in the sky are not so clear since many modes are competing for 
attention. The superposition of many otherwise distinct geometric patterns can lead to 
something apparently random. The task of geometric methods is to separate the patterns 
out of the sky. The statistic of eqn. ( |6.11C| ) manages to draw out circles. Other patterns in 
addition to the circles can emerge as described in this section. 



The patterns are best siphoned off a CMB map by scanning for correlations [15|. As an 
example consider a map of the antipodal correlation 

A{h) = CM{n,n) (6.119) 

which measures the correlation of fluctuations received from opposite points on the SLS 



|15|. In a simply connected space opposite sides of the SLS should have no communication 
between them and a map of antipody would generate nothing more than a monopole. It 
is possible that accidental correlations appear at random but no geometric structure would 
emerge. By contrast, if the universe is topologically connected, then points on the manifold 
which seem to be far apart in the tiling picture may actually be quite close together in the 
fundamental domain. Therefore opposite points on the SLS may be strongly correlated, 
may in fact be the same point. This is another example of ghost images but in an antipodal 
map collections of ghosts are caught and a picture of the symmetries of the space emerges 
PI. 



78 



Figure 43: 



The size of a spot can be estimated at the Silk damping scale below which fluctuations 
have smoothed some. The angular size of these spots are too small for COBE— DMRto have 
detected but will be visible to the high resolution MAP and Planck Surveyor. 

Although the search for pattern formation in correlated maps is model independent, a 
zoo illustrating the variety of structures compact manifolds produce can be built with some 
simple approximations. The correlation between two points on a compact manifold can be 
estimated as the correlation they would have on the universal cover given their minimum 
separation: 

CM(.n,h')^C^[d^in(.x{n),x'{h'))] , (6.120) 

where is the correlation function on the universal cover and dj^m is the minimum distance 
between the two points in the topological space. The estimate is effectively the lowest order 
term in the method of images. It is inadequate for use in a likelihood analysis but is sufficient 
for predicting the types of patterns which emerge from topological lensing. The images of 
a given point out to order m are found with the generators of the identifications as 

m 

= n5fc»^(^') • (6.121) 

i 

The image point which lands closest to x{n) determines dmin- 

As an example, we show the antipodal map for a 27r/3- twisted hexagonal prism in fig. 



44 For the antipodal map we prefer the orthographic projection which shows the genuine 
shape of the surface of last scattering instead of the Aitoff projection customary in 5T{n) /T 
maps. Notice the clear hexagonal face drawn out by the correlated map. Since antipody 
is symmetric under vr, the back is a copy of the front. For the vr/S- twisted hexagon, if the 
space is small enough, pairs of circles appear in A{fi)^ as shown in the right-most panel of 
fig. 44. These are the circles in the sky of section § |6.2.1| . 

A compact hyperbolic space shows distinct patterns. The compact icosahedron known 
as the Best space after the mathematician who identified the manifold provides the best 
testing ground |123|] . The map of A{n) is shown in the left of fig. 46. Antipody outlines 



pairs of identified triangular faces and also locates circles. Clearly a symmetry group for 
the Best space is located in this map. Another correlation function is also shown which 
compares one point on the SLS to the rest of the sphere. The point selected is a copy 
of the origin and so reflects the most symmetric observation of the fundamental domain. 



Another example is given by the antipodal map for the Thurston space in fig. 47. Some of 



the correlated features such as the arcs in fig. 47 may be secondary correlations and it is 
not clear they will ever be bright enough to be observed. 

In fairness, it is difficult to know if any of these patterns will really be measurable 
in a realistic experiment with physical complications such as the thickness of the surface 
of last scatter, additional Doppler effects, noise, fictitious correlations etc.. To read the 
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Figure 44: Left: Orthographic projection of A(n) for a hexagonal prism with L = O.6A77. 
Right: A(h) for 7r/3- twisted hexagonal prism. The length of the prism direction is .24 while 
L = 1. There are circles 




Figure 45: Dirichlet domain for the Best space. 




Figure 46: Left: A{h) for the Best space with Q = 0.3. Right: The point-to-sphere correla- 
tion. 
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Figure 47: Antipody in the Thurston manifold. 



correlations from the future data, real space statistics will need to be developed which 
handle smoothings, subtractions of low order multipoles, noise and fictitious correlations. 
While some statistics have been promoted, a realistic approach will likely develop only when 
the data is actually available. 
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Figure 48: An observer O at rest with respect to the compact spacehke hypersurface S 
and an observer O moving at constant velocity with respect to S on a periodic orbit of the 
torus. The orbit corresponds to Xend = TyT'^Xst&rt- 



7 Beyond Standard Cosmology 
7.1 The twin paradox and compact time 

So far we have ignored the issue of time. We have exphcitly considered spacehke hyper- 
surfaces S of constant curvature fohated by a natural conformal time. Compactification of 
these surfaces leads to the pictures we have described without compactifying time. However, 
an observer moving on an inertial worldline which is not at rest with respect to the cosmic 
expansion will perceive an identification which mixes spacetime coordinates as dictated by 
the Lorentz transformations A so that x = Ax with x = (rj, x) comoving coordinates. As a 
result, a compactification of S in the comoving coordinates of the form (ry, x) — > (r/, x + L) 
for instance will result is an identification of a time shift as measured by the non-comoving 
observer of (1 — f3'^)~^/'^ (ji — P ■ x^ ^ {1- /5^)^^/^ (jj — (3 ■ {x + L)^ where /? is the velocity 
relative to S. 

As an explicit demonstration, consider the twin paradox [124, 125 , 126 1. Let O be at rest 
in a comoving frame where a flat spacelike hypersurface is compactified into a torus. Let O 
be an inertial observer on a periodic orbit with respect to O as in fig. The periodic orbit 
will obey the boundary condition x = Ax — > A7X where 7 is the corresponding word. From 



2.3 we have 7 = TyT^ (to include the physical time in the embedded coordinates and in 



Tx, Ty see Ref. |124|] ). According to O, both space and time points have been identified. As a 



result it becomes impossible for O to synchronize her clocks |127|. The lack of synchronicity 



will be given by the time component of A(l — 7)x [124]. 

Since both O and O are inertial, by the principles of relativity, each should believe 
the other's clocks run slower and therefore each could expect the other to be younger 
at their reunion leading to a paradox. However, the topological identification breaks the 
general invariance and selects a preferred frame, namely the frame in which the topological 
identification is purely on S. In that frame clocks can be synchronized and the volume 
of space looks smallest. The observer on the periodic orbit, though inertial, will discover 
that their clocks cannot be synchronized and this additional shift leads them to compute 
an older age for their twin. Both agree the twin at rest with respect to S is older. 

This simple thought experiment emphasizes that compact topology selects a preferred 
frame, namely the frame in which the universe looks smallest and in which observers can 
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synchronize their clocks |127] (see also Refs. [128]). To generalize to curved space, A can 



be replaced by an appropriate diffeomorphism and the spacetime topology generalizes to 
Mc = R^G/T. 

Finite spaces reverse some of Copernicus' philosophical advances by selecting a preferred 
location at the center of the space, a preferred observer at rest with respect to the com- 
pactification and a preferred time. While Copernicus may have removed us from the center 
of the universe, topology puts some observers back there. 

This also raises the question of compactifying time outright. Probably time could be 
compactified so that there were closed timelike curves that were not causality violating. Very 
restrictive possibilities would result since only events which could repeat ad infinitum would 
obey the boundary conditions. It would be hard to envision a compact time model being 
consistent with the laws of thermodynamics, except perhaps on a cosmological timescale. 
The compact time scale would have to be much shorter or much larger than biological 
timescales or no children could sensibly be born since they would somehow have to grow 
young again. A universe could go through a big bang and eventual big crunch only to repeat 
the history of the cosmos with another big bang. The same galaxies, stars and the same 
people would be born, live and die. Fated to repeat their paths ad infinitum. Quantum 
gravity may reset the initial conditions at each big bang allowing new galaxies, new stars 
and new organisms to form, sparing us from a relentless cosmic boredom. The fanciful 
possibility of a compact time magnifies the already strange and distinct nature of time. 

7.2 Extra dimensions 

Topology it has been suggested is a discrete feature of space As such, it may be better 
integrated in a quantized theory of gravity while there is no prediction for topology in clas- 
sical relativity. The earliest attempts at creating a finite universe grew out of semiclassical 



quantum cosmology |12£, 130, 131, 132]. It may seem intuitive that a smaller universe 
should be easier to create from nothing than a larger universe. Therefore the probability 
for creating a small, finite cosmos may be relatively high, if only we could compute the wave 
function. However, technical and conceptual difficulties dog the semiclassical approach such 
as defining a measure on the space of states, normalizing the wave function etc. A convinc- 
ing statement about the topology of the universe will very likely require a fully quantized 
theory of gravity. 

Interestingly, additional dimensions have featured prominently in attempts to quantize 
gravity. These extra dimensions are always topologically compact and very small. While 
no quantum gravity theory is yet able to predict the topology of space, the possibility 
that compact internal dimensions will have topologically compact external (that is, large) 
dimensions is certainly alluring. Ultimately, a fundamental theory should predict the global 
topology of the entire manifold whether it be 3-D or IID. In the meantime the hierarchy 
between small and large dimensions remains mysterious although some recent suggestions 



have created a bit of a stir 133] 



Kaluza-Klein theories introduced extra compact dimensions in an attempt to unify fun- 



damental theories of physics ]134]. Upon compactification, the radii of the small dimensions 
behave as scalar and tensor field theories. More modern string theories naturally invoke 
extra dimensions in a manner reminiscent of these early Kaluza-Klein models. Recent fervor 
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in string phenomenology has involved compact extra dimensions of moderate to large size 
in an attempt to explain the hierarchy problem in standard particle physics. The hierarchy 
problem questions the disparity in scales from the Planck mass of 10^^ GeV to the mass 
of the electron at a few eV. A unification of fundamental theories has to naturally justify 



this span over 25 decades of energy scales. As in Ref. [ 135 | we consider a spacetime with 
4 + dimensions of the form Ai = R<Si (8) Ai''^ where R represents time, is a con- 
stant curvature Friedman-Robertson- Walker metric and Ai-'^ = G-^ /T is an AA-dimensional 
compact internal space. The (4 + A^) -dimensional gravitational action becomes 



A = ( S^+^^x,J-g(^+^)n^^+^h'll (7.1) 



In the simplest Kaluza-Klein picture, integration over the compact extra dimensions leads 
to the (3 + l)-dimensional action we experience 



A ~ / d^xJ-g(^)n^^^M^i + ... (7.2) 



plus additional dynamical terms where M^i = M'^^-'^ and R is the radius of the internal 
dimensions. The hierarchy between Mpi and is large if RM^ is large |135| ]. The disparity 
in energy scales then becomes a dynamical and geometric question. 

The modern ideas inspired by string theory involve the localization of matter to a 3- 
brane nested in the full (4 + AA)-dimensional space. The scale is expected to be ~TeV 
as predicted by supersymmetry. Since ordinary matter is confined to a 3-brane, we would 
be unaware of these extra dimensions regardless of their size unless we try very hard to look 
for them. Some laboratory experiments are underway to probe any additional dimensions 
which may be lurking there. 

The topology of the compact extra dimensions is not well understood. There are only 
a few requirements the internal dimensions must satisfy in order to break supersymmetry 
at a physically sensible scale. The internal spaces fall loosely under the broad category 
of Calabi-Yau manifolds. Nearly all topologies investigated so far are modeled on flat 
geometries. An outgrowth of cosmic topology has been the suggestion that these internal 



dimensions be compact and negatively curved |135|. The hyperbolic internal spaces have 
some advantages over flat space including a less demanding tuning of geometric parameters 
and a suppression of astrophysically harmful graviton modes from the extra dimensions 



1 135]. Another advantage of the compact hyperbolic extra dimensions over compact flat 
extra dimensions may be that chaotic mixing on finite hyperbolic manifolds could explain 
the smoothness and flatness of the large dimensions as argued in Ref. |^]. We briefly 
discuss chaos in the next section. 

The profound connection between small dimensions and large have only begun to be 
forged. As many times in the past, cosmology provides a unique terrain in which to test 
fundamental theories [ |136| | . 

7.3 Chaos 

We have touched upon the chaotic motions of particles on a compact hyperbolic space. 
Chaos refers to the thorough mixing of orbits which show an extreme sensitivity to initial 
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conditions. Geodesies deviate on a surface of negative curvature as mentioned in §^ so 
that they are extremely sensitive to initial conditions. The motion becomes fully ergodic 
upon compactification of the surface. Formally, chaos in these Hamiltonian systems means 
the geodesies equations are nonintegrable; there is no smooth analytic function which can 
interpolate between orbits with different initial conditions. 

Despite the resistance of chaotic systems to conventional integration methods a great 
deal about the structure of phase space can be determined. Much like thermodynamics, 
chaotic systems can be understood in terms of a set of states dense in the phase space. 
For chaos this set is provided by the collection of unstable periodic orbits which grow 
exponentially with length. All aperiodic orbits can be understood in terms of this special 
subset. The periodic orbits pack themselves into a fractal set in order to fit within the finite 
phase space. Recall the fractional dimension maintains the set at a finite volume but allows 
for an infinite area. In this way fractals try to maximize the information content, so to 
speak, while minimizing the volume. The explicit fractal structure on a compact octagon 
was isolated in Ref. [124] 



Pursuing the analogy with thermodynamics, entropy can be defined as 

N{k) 

M/x)= lim - ^ /xln(l//.(A;)) (7.1) 

K — y OO rC 

I 

The calculation goes like this. Draw a fundamental domain. Now use the generators gi.-.gn 
to tile H'^ with copies. The measure at order k is then defined as the fraction of the total 
volume where each tile has the same volume: 

= ^ (7^2) 
and N{k) is the number of unique tiles. The sum in eqn. ( [7. ID reduces to 

h{pi) = lim - \n{N{k)) = Ht- (7.3) 

k^oo k 

This is equivalently the topological entropy, a symbolic entropy which counts the number 
of accessible states; that is, closed loops or equivalently tiles in the tessellation. 

Recall that periodic orbits can then be counted symbolically with the homotopy group. 
We do not distinguish between loops of varying lengths if they are homotopic. The spectrum 
of periodic orbits so defined are therefore a topological feature. The number of unique words 
that can be built of length k out of a n-letter alphabet with r relations is equivalent to the 
number of neighbors at order k in the tiling. Not all of these neighbors will be unique. The 
repeats are accounted for by the relations and we know that the number of neighbors at 
order k is bounded by 

{2n)k > N{k) < (2n)*^ (7.4) 

the lower bound being simple spaces such as the flat nonchaotic torus and the upper bound 
the unpruned maximally chaotic case. The topological entropy of the torus is zero while 
the entropy of an unpruned hyperbolic space is Ht = ln(2n). The compact octagon was 



shown to have a topological entropy of Ht = In 7 |124| 
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Another kind of entropy, the metric entropy introduced in describes how quickly 
mixing takes place in a chaotic system and can be expressed as the sum of the positive 
Lyapunov exponents: h = J2'^- 

This is just a taste of a rich area in dynamical systems theory and the reader is referred 
to the many excellent texts on the subject |90| . 

Although cosmologists gingerly avoid this complex feature of the dynamics, chaos does 
have profound and unavoidable consequences. One of the original motivations for a compact 
topology exploited the chaotic motions. The chaotic mixing could lead to a dilution of initial 
anisotropies leading to the symmetric universe we observe today. These initial attempts 
failed since the space could not be made small enough to allow for sufficient mixing and still 
be the large cosmos we observe today . Variants on this idea fuse the chaotic mixing with 
an J7 < 1 inflation model where an early episode of chaotic mixing provides the moderate 
initial conditions needed to permit a subsequent inflationary phase to succeed |4^. 

We have already encountered the nonintegrability in trying to find eigenmodes on com- 
pact hyperbolic cosmologies. For the initial spectrum of fluctuations, researchers continue 
to use the assumption of a flat Gaussian spectrum as motivated by inflation. Still, since 
the consistency of inflation and an observably small topology is shaky, a more consistent 
approach would be to examine a distribution of initial fluctuations on a compact space 
in the absence of inflation. The numerical results of §6.1 do suggest that a Gaussian flat 
spectrum is in fact natural on compact hyperbolic spaces, even in the absence of inflation. 
While this work is very suggestive, a primordial quantum system has never been very thor- 
oughly thought through. Because of the importance of this issue we take a moment to 
discuss quantum chaos. The quantum system is relevant to our discussion regardless since 
the expansion of fluctuations will be analogous to the expansion of the semiclassical wave 
function. The stationary Schrodinger equation is the usual Helmholtz eqn ( |6.78| ). 

The quantization of the chaotic system is still not fully understood but interesting 
features have been conjectured and confirmed. Reminiscent of the Feynman path integral 
approach, the wave function can be thought of as the sum of classical trajectories. Since the 
classical trajectories chaotically mix in phase space, it has been conjectured that the wave 
function would be a random Gaussian function of the eigenvalue q with a spectrum given by 
the Wigner distribution function. Wigner's function is an attempt to generalize Boltzmann's 
formula for the classical distribution of a statistical system to a quantum system [137|. It 
provides a description of general attributes of wave functions 1 90 1 and shows how the function 
tends to distribute itself over classical regions of phase space. Notice the similarity between 



this suggestion and the method-of-images as well as the numerical results of Ref [ [100| , 103 

The nature of the discrete spectrum of energy levels can also be related to the classical 
chaos. As already argued in the numerical work on the eigenmodes of a compact space, the 
number of eigenstates can be related to the volume of the space through Weyl's asymptotic 
formula. Additional features, such as the spacing between energy levels, can be directly 
associated with the underlying chaos. 

Despite this conjecture of randomness, distinct remnants of the underlying chaotic dy- 
namics have been found in the form of scars [138|. Scars are regions of enhanced intensity 
in high k states along periodic orbits. It is hard to imagine the very high k modes having 
cosmological significance since it is the low k modes which probe the largest cosmological 
distances. Scaring is expected to be less prominent in low k modes since the width of the 



86 



enhancement will be correspondingly more diffuse. Although diffuse, scars on large scales 
could provide the small catalyst needed to order structure on large-scale. Scarring on a 
2D double doughnut was investigated in references |88|. The underlying fractal structure 
of periodic orbits on the double doughnut was studied in detail in Ref. [ [L24| ]. There it was 
suggested that the filamentary structure we observe in the distribution of galaxies may be 
a consequence of this slight enhancement of the seeds of structure formation. 
For other important chaos articles see Refs. 142 , 141, 140, 139 1 
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8 Summary 



The creation of the universe is still not well understood. Clarity on the earliest moments 
will likely come only with a fully functional quantum theory of gravity. In the meantime, 
we know that space is curved and evolving and most also possess some topology. If the 
curvature of the universe falls within the observable horizon, then topology may also. The 
CMB provides the deepest probe of the universe on the largest-scales and we have reviewed 
the many ideas on how to extract the topology of space from maps of the microwave sky. The 
methods fall into two primary categories: direct statistical methods and generic geometric 
methods. The salient features in CMB maps which reveal topology are (1) a discretization 
of the sizes of hot and cold spots, (2) a cutoff in the spectrum for wavelengths too big 
to fit within the finite space and (3) an anisotropic and inhomogeneous distribution of 
correlations corresponding to repeated ghost images of the same spots. Data from the 
future satellite missions MAP and Planck Surveyor are needed to determine if the CMB 
docs in fact encode such features. When the new high resolution maps are in hand, we 
will have to face the potentially prohibitive difficulty of foreground contaminations, the 
thickness of the surface last scattering, and fortuitous correlations. If we are lucky enough 
to surmount these observational trials, we may be able to see the entire shape of space. 

Still, as many fear, the topology scale may naturally be far bcj/'ond the observable 
universe. If this is the case, we can turn to physics on the smallest scales to learn something 
about what we will never see on the largest scales. If an ultimate theory of gravity beyond 
Einstein's is able to predict the geometry and topology of small extra dimensions, there is 
every reason to hope we will learn, if only indirectly, the geometry and topology of the large 
dimensions. 
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A Appendix: Representations of Hyperbolic space 



Because of the recent emphasis on hyperbohc models, we include a section here on other use- 
ful coordinate models for manipulating H^. We have already seen two different coordinate 
systems for the metric eqns. (|4.2| ) and ( [4.4[) . Another important coordinate system comes 
from embedding in the (3 + l)-Minkowski space as discusssed in §2.1. The Minkowski 
metric is 

da'^ = -dxl + dxl + dxl + dx| (A.5) 
The 3Z)-hypersurface is constrained to have pseudoradius 



Xq + + X2 ~\~ 



1 



(A.6) 



and the curvature radius has been taken to 1. In Minkowski coordinates the generators 
7 € r take the convenient form of 0(4), orthogonal 4x4 matrices. These coordinates can 
be related to those of eqn. (4^) with the transformation 



x"- 



/x^\ 

^1 



\x') 



( 



cosh X 
sinh X sin 9 cos ( 
sinh X sin sin <: 
sinh X cos 9 



(A.7) 



There are several other useful representations of the hyperbolic plane including the 
Poincare ball with 

X = tanh(r/2)n (A.8) 

with h the usual unit vector in spherical coordinates. The Poincare ball model maps H'^ to 
the open ball {x S E^|x • if < 1} with 

4dx • dx. 



da' 



1 



X • XI 



It is useful to know the geodesic distance is 



d{i, x') = arccosh 



1 + 



2|x- 



X 



(l-|x|2)(l-|x'|2) 



(A.9) 



(A.IO) 



All geodesies intersect the boundary orthogonally and are therefore semicircles or straight 
lines which are the diameters of the ball. 

There is also the upper-half space representation {x G E'^jxs > 0}, 

(dx^ + dy^dz^) 



da' 



z- 

dp" + e-^P [dx^ + dy-") 



(A.ll) 



with e'' = z. This coordinate system is particularly useful for cusped manifolds which are 
discussed in § 6.1.2| . 

Finally there is the three-dimensional Klein model with coordinates 



tanh xn 



(A.12) 



which is used often in constructing the Dirichlet domain in three-dimensions and mapping 
the periodic geodesies. Geodesies are mapped into straight lines in this model. 
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